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CHAPTM I 

INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

The Augmentor Wing Jet STOL Research Aircraft CAWJSTOLRA) 
is an advanced test vehicle for STOL powered- lift flight 
testing. An initial area of interest is the control of 
the AWJSTOLRA to closely track an altitude profile during approach 
and flare maneuvers. For a powered lift vehicle during landing 
approach, the throttle may be used as the primary control of 
normal acceleration. The accurate digital simulation of 
aircraft maneuvers involving throttle changes and the design 
of autopilot logic using feedback to the throttle require an 
accurate model of the dynamic response of engine thrust to 
changes in throttle commands. This model represents a complex 
interaction of nonlinear static and dynamic behavior in the 
engine control and aero thermodynamic and mechanical couplings 
in the engine and duct. It has been th^e goal of Systems Control, 
Inc. 'CVt) [SCI CVt)], in close coordination- with NASA Ames personnel, 
to determin'e experimental mathematical models of the AWJSTORA 
Spey engine dynamics by the application of system identification 
techniques to flight test data. 

The AWSTOLRA is a modified de Havilland Buffalo C8-A air- 
craft using two Rolls Royce- Spey engines [1]. The program is 
sponsored jointly by the National Aeronautics and Space Adminis- 
tration and the Department of Industry, Trade and Commerce of 
Canada. 

Modifications to this C8-A were made to improve lift 
characteristics at low speed for the short landing requirement. 

A transverse duct is used to channel fan airflow to ejector 
nozzles on the opposite wing blowing over an augmentor flap. 

Direct lift is provided by vectoring core stream flow 
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through an adjustable nozzle. The modifications have a signi:^ 
ficant effect on the dynamical properties of the propulsion 
system and lead to the consideration of experimental model 
development for this aircraft. 

The Spey engine is a low bypass, twin spool turbofan built 
by Rolls Royce of Canada. The normal gas path has been signifi- 
cantly modified for AWJSTOLRA operation. The new configuration 
is shown in Figure 1.1. The engine fuel control is hydroraechanical 
with a mechanical cable linkage used to command throttle lever 
inputs from the cockpit electromechanical servo. 

The linkage/control/engine system is nonlinear in both 
static and dynamic operation. Simple models previously utilized 
to describe engine behavior were not sufficient to explain ob- 
served response to closed- loop throttle inputs during landing. 

The SCI (Vt) program was undertaken to utilize system identifi- 
cation software to develop flight verified models of propulsion 
system dynamics. 

1 . 2 SUMMARY 

The report contains six chapters describing the develop- 
ment and validation of the Spey engine model. Chapter. II de- 
scribes the analysis of the dynamical interactions involved in 
the propulsion unit and the reduction of the model to contain 
only significant effects. The reduced model was used in con- 
junction with the flight data to develop initial estimates of 
parameters in the system using classical approaches. This 
analysis is presented in Chapter III. Chapter. IV describes 
the theoretical background used in estimating the parameters 
and the -software package developed to accomplish the processing. 
Flight data was processed to estimate and validate a set of 
parameters and these results are presented in Chapter V. The - 
v;ork and results are summarized in Chapter VI. 
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Figure 1.1 Gas Path Configuration of Spey Engines on AWJSTOLRA 


CHAPTER II 

MODELING. THE SPEY ENGINE 


2.1 INTRODUCTION 

A simulation of th.e AWJSTOLRA was required to assess its 
handling and control characteristics prior to flight test [1,2]. 
In addition, pilots were to receive training in various flight 
regimes including system failures and abnormal operation. For 
this reason, a detailed digital simulation of the AWJSTOLRA 
aerodynamic and kinematics was developed by NASA Ames, including 
a model of the characteristics of the propulsion system. In 
this chapter, models for the propulsion system of varying levels 
of complexity are examined and a model is proposed for system 
identification . 

2.2 TYPES OF ENGINE MODELS 

The propulsion system can be described in mathematical terms 
in several ways. The choice of model form has an important im- 
pact on the analysis tools, the type of data taken, and the 
regimes of application of the expressions. Another consider- 
ation in the model is the interface between the engine and the 
input /output boundaries. Models that describe the throttle 
linkage and fuel control as well as the engine may be signifi- 
cantly different from "open-loop" models of the engine . 

Also, the estimation capabilities from experimental data will 
be directly influenced by the model interface. This partition 
of the physical system into model and "outside" is determined 
by the objectives of the descriptions. • 

In the Spey engine, a model is- desired which describes the 
engine thrust and airflow response to power lever inputs from 
the cockpit. Dynamics in the engine or control which are not 
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excited during tliis operation are not important to the modeling 
task. The model will be used to describe large power excur- 
sions approximately and model dynamics in the high power 
region accurately. Altitude/speed regions o£ validity will be 
limited to low altitude/low speed and nonstandard temperature 
conditions . 

The value o£ the model will be in its description of 
important aspects of performance. Steady- state values of 
thrust and airflow should match the values measured in static 
tests corrected to nonstandard conditions. Dynamic response 
to small perturbations should accurately reflect engine behav- 
ior. ■ Match of intermediate variables describing engine 
operation is not critical , 

Models used to describe typical engine behavior vary 
widely in complexity and accuracy . [3] . The simplest descrip- 
tion of engine response is to plot corrected engine thrust 
versus power lever angle. This characteristic is shown in 
Figure 2.1. The'plot can be used as a dynamic model if an 
appropriate time lag is associated between the actual throttle 
position and a lagged or virtual position. The engine is 
observed to accelerate and decelerate at rates dependent on 
the power level. This effect can be represented in the model 
as a variable rate limit. The experimental model is shown 
in Figure 2.2. No attempt is made at a phenomenological 
explanation of the behavior and parameters are adjusted from 
observations. Since the model does not reflect the internal 
derivative agents of the thrust response, the match between 
the system behavior for various size inputs and for different 
starting and ending values of thrust will be poor. A more 
significant drawback of these models is the poor closed- loop 
description of the behavior. This mismatch is typically 
manifested by an incorrect prediction of the closed-loop 
stability of the system. This is due primarily to matching the 
step response in estimating parameters rather than matching 
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Figure 2.1 Thrust Response to Throttle Input 





Figure 2.2 Experimental Model o£ Engine Dynamics 











frequency response over a suitable bandwidth. Generally, 
frequency matching requires more complex model forms. 

In order to accommodate this lack of precision, simple models 
tend to be ’’improved*’ by ad hoc additions which attempt to 
correct the fundamental nonphehomeno logical characteristics 
of the formulation. 

Far more complex digital simulations exist which can 
model component characteristics measured from rig tests and 
aero thermodynamic phenomena occurring in the gas path. Such 
a simulation is shown schematically in Figure 2.3. These 
simulations include basic physical laws relating energy, 
work, massflow and acceleration as they dynamically interact 
in the engine. Various ’’adjustable" parameters such as lumped 
isentropic efficiencies and areas are adjusted to match the 
observed steady- state relationships between input and output 
variables, A detailed representation of the governor is 
included and the overall response can be generally tuned to 
match observed dynamics.; ‘TVhile this type of simulation 
recreates steady-state performance accurately, the match of 
transient response is poorer due to modeling uncertainties 
in the complex equations. Predicted stability characteristics 
in a closed- loop control are quite good. Thus, these programs 
can be used as a durable test bed for control and autothrottle 
design and evaluation. 

The drawbacks of detailed analytical simulations are 
that: 3- large computer capability is necessary to execute 

them (often slower than real time) , (b) detailed internal 
engine and control information is required which is only 
available Gif at all) to the engine manufacturer, and (c) the 
added complexity cannot be justified or validated from input/ 
output performance observed in operation. Clearly, a middle 
ground must exist. 

The development of an experimentally validated analytical 
description of the propulsion system can be accomplished by 
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considering the available measurement data and the type of 
inputs provided, Tke development starts with a complete analysis 
of possible dynamic effects. Data is examined and unimportant 
or unexcited phenomena are removed. The remaining model terms 
are examined for bandwidth of response relative to the simu- 
lation requirement. High frequency effects can be neglected. 

The resulting system represents a description of the response 
including both steady-state and dynamic performance whose 
parameterization is available from recorded measurements. 

This model will provide an excellent component module in an 
overall aircraft simulation used as a preliminary tool in the 
design of outer loop integrated control functions such as 
flight direction, autothrottle and SAS. 

2.3 CURRENT SPEY ENGINE MODEL IN THE AWJSTOLRA SIMULATION 

The engine model used on the AWJSTOLRA simulator is shown 
in Figure 2.4 C^- detailed listing is included in Appendix B] . 

The throttle setting at the input to engine control, 
differs from the commanded value by a hysteresis function 
shown^ which models the mechanical backlash. Two forward 
paths model accelerations and decelerations. This form is used 
because the engine control accelerates the turbine using a 
different control law than the one used during decelerations 
as is discussed in Section 2.4. The forward paths modulate 
variable rate limits for the first order lag system which 
produces the lagged or virtual throttle position, ‘^VTH’ 
the output. The acceleration rate limit is low at low power 
and increases at higher power levels. This increased rate 
limit is lagged by the fuel metering time constant, t^, which 
is discussed in Section 2.4. The deceleration path has a 
small signal time constant which is 60% of the acceleration 
time constant. This represents a parameter used to match observed 
full throttle transients. It is not internal to the engine, 
since for small inputs, the engine acceleration and deceleration 
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Figure 2^4 



Spey Engine Model in AWJSTOLR Aircraft Simulation 












time constants are the same, Stated another way, th,e engine 
can be modeled as a l-inear, dynamic system for reasonably 
small inputs. The deceleration rate limit decreases at 
lower powers. Figures 2.5 and 2,6 show typical sea level 
static acceleration/deceleration time responses. Note that 
the form of the model in Figure 2.4 is dictated by the response 
rather than fundamental system phenomenology describing the be- 
havior. 

It can be observed from the model that, at power levels 
above about 65% thrust, the engine is operating away from the 
rate limits for most inputs. Removing these blocks from the 
figure produces a simpler view of the equivalent high power 
representation of the system. This is shown in Figure 2.7. 

The system is quite simply represented, but is fundamentally 
nonlinear since the loop gain and time constant depend on the 
sign of the power’ transition. 

The model in Figure 2.7 uses steady-state representations 
to calculate the mass £low_ and thrust. This is equivalent t'o. 
the assumption that the engine maintains its steady-state 
rotor speed match during nonsteady operation. This assumption 
is not correct. Cold and hot thrust are determined by the fan 
airflow and core airflow which have independent degrees of 
freedom. 

In summary, the possible problems in utilizing a model 
of the form of Figures 2.4 and 2.7 are as follows: 

(1) A nonphenomenological representation of behavior 
leads to false conclusions such as a fundamental 
nonlinearity of response. 

(2) A small number of degrees of freedom produce poor 
transient matches. 

C3) Matches of engine response between arbitrary end- 
points are uncertain and the number of parameters 
in the model is insufficient for adjustment. 
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Figure 2.5 Engine Test Stand Acceleration 

Characteristics 
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Figure 2.6 Engine Test Stand, and Flight Test 
Decelerations to Idle 
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Figure 2.7 Simplified High Power Model 








2.4 DETAILED MODEL OF THE ENGINE 


Tlie Spey engine is characterized by a group of nonlinear 
dynamical elements; Ca) the engine, Ci») the fuel control, and 
Cc) the throttle linkage. Models representing each of these 
elements can be derived for the individual physical members 
in the system. Simple analysis of each component produces 
approximate parameter values. The model can be simplified 
for identification to contain only those elements which can 
be estimated from the flight data. The development of the 
detailed engine model will be described in this section. 

The dynamics of the engine in the region near full power 
are nearly linear. Small perturbations in the state of the 
engine are characterized by the two rotor speeds. Previous 
analysis of turbofan engine dynamics has verified the linearity 
and the structure of the linear system representing this be- 
havior [3,4]. 

Since the engine has no variable geometry in the gas path , 
the primary, dynamic exciter is the fuel' flow which is metered 
by the hydromechanical governor. Fuel is chemically converted 
to heat in the combustor. This provides excess energy at the 
turbine entrance. The combustor lag is typically on the order 
of 0.5 to 2.0 msec and is therefore much faster than the 
dominant system response. The gas temperature at the entrance 
to the turbines is converted into work by the expansion through 
the rotating turbines. This energy is used to drive the 
compressor on the high rotor shaft and the fan on the low 
rotor shaft. The dominant engine dynamical states can be 
associated with the rotating inertia of these elements. 

Newton’s law for the shafts can be written in terms of the 
speeds and torques as follows: 
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^^HTURB ' ^COMP^ 


C2.1) 


^^LTURB ~ ^FAN^ 

where C C refer to the high pressure and low pressure com- 
ponent shafts. ,The torques are determined from the mass flow and 
pressure change, across the component. The form of the expression 
is given below (for a compressor) 

"^in^CPr^'^/^-D 

Q ^ (^2.3) 

c 

where the compression ratio, is determined by flow equilibrium 

through the engine and duct. 

These equations can be used as the basis of a global nonlinear 
model of the engine. Nonlinear maps and detailed dynamical repre- 
sentation of nonequilibrium flow can model the continuity, energy 
and mixing phenomena occurring in the gas path. These terms 
can be considered negligible if one .is iinterested in system 
bandwidths of flight control, engine and autothrottle func- 
tion. 

For modeling the engine behavior during moderate speed 
excursions, a linear set of dynamical equations provides an 
accurate representation. The form of the system can also be 
simplified to second order for most turbofans by considering 
the physical character of the motion. For perturbation of 
fuel flow with bandlimited content, the two engine spools will 
accelerate in a collective fashion. For perturbation in fan 
loading due to variations in duct pressure or inlet distor- 
tion, the fan speed will tend to "relax" toward equilibrium in- 
dependently of the compressor. Using these observations, or 
equivalently calculating the relative magnitudes of the com- 
ponent map slope terms [5] such as 

IS 

3N » 

the qualitative analysis can be verified. 
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The linear system of constant ■ coefficient equations which 
is constructed for the engine is shown in a Laplace transform 
block diagram in Figure 2.8 and collected in state variable 
form as follows: 



Figure 2.8 Linear Engine Model 



- ° - 



“^FH^L 
- ^FH 


6W^ (2.4) 


where the five linear parameters, x^^, ^FH’ 

represent the dynamical behavior of the engine. 


and 


a 


The dynamical equations are written for perturbations 
away from an equilibrium condition. A typical choice for this 
equilibrium point is 100% power. 

In the Spey configuration on the AWJSTOLR aircraft, the 
ducting arrangement can cause the engine dynamics to be more 
uncoupled than the normal mixed flow path configuration. The 
fan flow is completely ducted to the active lift devices on 


18 







the wings. Back pressure effects on this flow will greatly 
affect fan response characteristics. The core stream exhaust 
is modified by a colander plate at the low pressure turbine 
exit. This plate provides choked exit conditions at all 
flight/power points and nozzle attitudes. The effect of nozzle 
motion on the flow is greatly reduced by this device. 

The flow equations modeling the pressure phenomena can 
be written from the continuity equations as follows: 

p = Am (2.5) 

where Am is the difference between inflow and outflow rates. 
This equation represents a flow lag effect and for small 
changes in power, the equation for the pressures in the duct 
and nozzle can be written: 

6Pd = CfiPp - 

^ tT “ ^EH^V 

£ 


( 2 . 6 ) 

(2.7) 


where the temperature dependence is ignored and the mass flows 
are modeled as linear functions of the perturbational rotor 
speeds which feed the volumes. Unlike the normal turbofan 
configuration, there is no flow coupling back through the duct 
to the fan from the nozzle volume. It is expected that the 
flow time constants will be small compared to the rotor lags, 
but because of the usual configuration, they will be retained 
in the model. 

For power levels typical of the approach configuration, 
the Spey engines themselves do not represent the largest 
source of nonlinearity. The throttle linkage and hydromech- 
anical fuel control represent the limiting elements in the 
linearity of the response. The components are quite complex 
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and accurate representations o£ the internal dynamics are dif- 
ficult to reconstruct from flight data. 

Engine power requests from the cockpit are transferred 
by a'^mechanical linkage driven by an electromechanical force 
servo. The cable moves the throttle arm at the engine. About 
a 15° rotation of the linkage represents the full idle - to 
maximum forward power excursion. The system contains a 
moderate amount of hysteresis which affects the small inputs 
typical of the autothrottle landing commands. Figure 2.9 shows 
the representation of the throttle linkage and the nonlinear 
equations are written below. This calculation is performed at 
each update time in the simulation. 



Figure 2.9 Throttle Linkage Model 


/ 6^(n) - h 5^(n) > + h 

6^(n+l) = < 5^(n)-h _< £ 5^(n)+h 

( 5-Cn) + h 6„(n)-h > 5 (n) 

^ C2.8) 

The hydromechanical fuel control is the most complex non- 
linear element of the system. Figure 2.10 shows the speed control 
and high pressure CH.P.) fuel pump. These systems affect the 
dominant response of the governor. A description of the 
elements of governor operation is included to put the detailed 
analysis into perspective. 

The fuel control meters fuel to engine in response to 
throttle inputs and engine demand. Instability, temperature, 
speed, and airflow constraints are accommodated while pro- 
ducing specified performance at different temperatures, 
altitudes and speeds. 
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Fuel is supplied from the tank via a low pressure pump, fil- 
ter, and flow meter to the high pressure CH.P.} pump entrance. Fuel 
is metered through the control to the fuel spray nozzles in the 
combustor, The pressure drop across the governor controls the 
fuel rate and this quantity is the primary actuated variable. 
Airflow limits, fan overspeed, temperature limits and emergency 
cutoff functions are performed by elements which are peripheral 
to the governor and not significant in the approach regime 
at near sea level conditions (the engine is not temperature 

limited at this condition 


The dynamical elements of the governor consist of the 
bellows assembly, speed governor, pressure drop governor and 
H.P. pump. Each element will be discussed and the appropriate 
model given. Figure 2,10 shows a schematic diagram of the 
pump and fuel control. 

The main element of the fuel system is the H.P. pump. 

This variable displacement device is constrained to maintain 
the pressure drop across the regulator by a servo piston feed- 
back.ielemeqt on the piamp stroke. The system is shown in Figure 
2.11. The linearized equation for this system is as follows: 

= 1/b [ACPy-P ) - kCW^-W^ 31 C2.9) 




A/b 



s + k/b 



Figure 2.11 Fuel Pump Model 


where the the pressure drop across the governor, 

and b, k, and describe damping and spring restraints 

on the piston. Fuel is metered by a pressure drop governor 
through a triangular orifice. The pressure drop across an ori- 
fice in a liquid flow is proportional to the square of the flow 
rate. By making the pressure drop proportional to the squares of 
the speed, a linear control governor action is achieved. A 
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Figure 2.10 Engine Fuel Control and Pump 
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Watt's governor can be used for this purpose with the rotating 
masses suspended- in the fuel to provide damping and fuel density 
compensation. The details of the mechanism are complex and a 
simplified diagram is used for the analysis. 

The triangular orifice is dra-wn in Figure 2.12. The metering 


ORIFICE 

SLEEVE 


Figure 2.12 Fuel Metering Orifice 



area is determined by two displacements, and X^. The 

X^ displacement measures motion of the throttle sleeve controlled 
by the governor and throttle springs. The position of the ori- 
fice is adjusted by; compression ratio and altitude by the cap- 
sule assembly discussed below. The metering area can be written 
as follows: 

\ C2.10) 

Considering only motion of the throttle sleeve, 

X^ = l/b(f^ - fg) (2.11) 

where f^ is determined by the spring action of the throttle 
linkage and f^ is the restoring force due to rotation of the 
counterweights. These may be xvritten as follows; 

£.j, = a^ + (X^ - X^) (2.12) 

fg ^ (2.13) 

where X^, is the commanded sleeve position derived from the 
throttle inputs. 
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The pressure drop across the metering orifice is sensed 
by the pressure drop pis-ton in the control. This piston moves 
an orifice in conjunction with a second set of fly weights. 

The orifice configuration is shown in Figure 2.13. 



Figure 2.13 Pressure Drop Governor Orifice 


and the equations of motion for this system can be written as 
follows: ■ 





. C2.14) 


and the total orifice area is given by the following 



and 


(2.15) 



C' 

Pd 


(2.16) 


The function of the governor is to provide stable, iso- 
chronous regulation at a wide range of operating conditions. 
Primary orifice size is limited because of dynamic range con- 
siderations. The operation can be understood by considering 
the static behavior. 

In steady state, the pressure drop governor has the 
following characteristics. ' 
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C2.17) 



Note that for fixed metering error (Eq. (2.19), = (2) the 

control law is unstable, i.e. increased speed increases the 
fuel. This occurs during' large transients-. when the governor 
area is limited by the throttle stops. Figure 2.14 shows a 
typical large -transient in the control plane. For_a_step -input, 
fuel flow increases at the pump response rate to the constant 
area line Fuel and speed increase along this line until 

the required speed is reached by the selected throttle curve. 

The throttle sleeve moves off the stop and the cubic speed 
relationship dominates the control law. The system is stable 
as shown. Accelerations follow the reverse logic. This method 
provides fast and stable fuel regulation for the Spey except 
at low power. Altitude variations are not compensated by the 
components. Note that antisymmetric rate limits will cause 
the engine to accelerate and decelerate at different rates for 
large, throttle limited power excursions. 

The capsule assembly and limiter provide altitude 

compensation and low power acceleration limiting for the engine. 

A model of this system which affects the orifice sleeve position, 
X 2 , and the metering area, is not described in detail. The 
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Figure 2.14 Governor Characteristic 

effect of these components is most pronounced during acceler- 
ations from low power. At these conditions, the fuel flow 
increase is slowed considerably until the rotor speeds are 
increased. This limits the possibility of surge during acceler- 
ation. Hysteresis can be important during small motions of the bel- 
lows and this might affect the overall repeatability of the response 

The various component models have been interconnected 
in Figure 2.15. The dynamic range of many of these elements is 
higher than the dominant response time constants. Thus, while 
the dynamical relationships can be ignored, the static feed- 
back structure of the governor will significantly ■ affect 
engine response times and linearities. 


26 




rigure 2.15 Detailed Dynamics Simulation of Spey Fuel Control 



















2.5 SIMPLIFIED MODEL- 


The system model can be simplified to include the linear 
model of the engine’, the hysteresis model of the throttle linkage, 
the variable metering orifice characteristic expanded to second 
order and the pressure drop regulator characteristic. This 
system is shown in Figure 2.16. 

The time constants of the electromechanical servo ar'e 
between 0.05 and 0.1 seconds which is faster than the ex- 
pected closed- loop engine response. Since commanded throttle 
and actual deflections are recorded, any relative time delay 
can be detected. The fuel pump response time should be faster 
than 0.1 seconds at high displacements. This response time 
is the limiting element in the overall governor response loop 
since fly weight motion should have response times less than 
25 msec. The engine control reduces to a nonlinear feedback 
regulator driving a rate servo modeling the action of the 
pressure drop governor. -Any mismatch between the forward and 
feedback path gains in the pressure drop system will result 
in a slightly stable or unstable servo loop. This gain 
mismatch can be lumped with the speed governor characteristic 
without loss of generality. The forward and feedback gains 
in this loop are assumed identical. The closed-loop character 
istic of the pressure drop governor is a pure integration (i.e. 
a rate servo on the speed governor output) . The nonlinearity 
of the metering orifice is reflected by a slightly nonsym- 
metric transfer characteristic. The overall hysteresis of 
the control after the throttle is represented by a second 
hysteresis block. 

A further simplification of the control model is possible 
assuming that the transfer nonlinearity and control hysteresis 
are small. This is discussed further in Chapter V. 
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2.6 SENSOR MODELS 


Data on Spey performance was acquired in flight. Instru- 
mentation was used as listed in Table 2.1. Speed, position, and 
pressure transducers were assumed to respond faster than the 
dynamics being estimated. The measurement errors associated 
with these inputs can be considered uncorrelated random pro- 
cesses with RMS values shown in the table. 


Table 2.1 
Sensor List 


VARIABLE 

SYMBOL 

RMS ERROR 

Low Rotor Speed 

\ 

0.5% 

High Rotor Speed 

"h 

0.5% 

Throttle Command 

5 

c 

0.2 Degrees 

Throttle Angle 

"^T 

1.0% 

Duct Pressure 

Pd 

2.0% 

Exhaust Pressure 

Pe 

2.0% 

Fuel Flow 

«fm 

2.0% 


The fuel flow was measured at the exit of the low pres- 
sure fuel pump. The sensor is a restrained flow turbine whose 
angular deflection is proportional to the turbine torque and 
mass flow. 

The system response is lightly damped because of the small 
value of flow damping supplied by the turbine. The model for 
the flowmeter is shown in Figure 2.17. 
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Figure 2,17 Linear Flow Meter Model 

The values ? = 0.2 and w =13.0 sec ^ were estimated from 

n 

the flight data. 

2 . 7 SUMMARY 

A group of mathematical models has been developed for 
the dynamical behavior of the Spey engine. Simple non- 
phenomenological forms are used to match specific transients. 
Detailed representations of the intercdmponent responses can 
be used to accurately predict the system response if dimen- 
sional and rig data is available. A simplified version of 
the detailed model is derived which contains important linear 
and nonlinear elements of the engine, control and throttle 
response for the landing approach power) configuration 

of the Spey. 
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CHAPTER III 

PRELIMINARY PARAMETER ESTIMATES 


The AWSTOLR aircraft was flown at Crow's Landing, Cali- • 
fornia to acquire engine response data for parameter estima- 
tion. The data was initially evaluated to determine initial 
parameter estimates in the dynamic model developed in Chapter 
II. Important terms were identified here and data consistency 
was validated. A nonlinear system identification routine was 
Incorporated with a digital- simulation of the engine. This 
program is described in Chapter IV. The results of the pro- 
cessing of the flight and the parameter estimates are compiled 
in Chapter V. 

3.1 SUMMARY OF FLIGHT DATA 

Flight data was taken in the approach configuration in a 
powered descent from 7500 ft. The 'nozzles in this configuration 
are rotated downward and the engines are operating near full 
power. The flight computer was used to generate throttle 
rate commands to the electromechanical servo. The engine 
throttles were modulated in unison. Appendix B contains a 
description of the flight log for the two data acquisition 
runs. Flight time for the two runs was 124' minutes from takeoff 
to final landing. Data was recorded during this time on 84 
instrumentation channels comprising engine, aircraft and avionics 
equipment sampled every 20 msec. Data records were chosen 
from the complete flight data for the parameter estimation. 

These data records are listed in Table 3.1. A set of instru- 
mentation inputs was chosen to reflect the variables incorpor- 
ated in the model. These measurements are listed in Table 3.2. 


PRECEDING PAGE BLANK NOT FMEn 
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Table 3.1 

Flight Data Records Used for Identification 


RECORD 

START TIME 
(d-h:in;s:ms) 

END TIME 
(d-h:m:s:nis) 

DURATION 

(SEC) 

NO. OF 
POINTS 

■BQH 

INPUT 

DESCRIPTION 

RAMP RATE 
(DEG/SEC) 

/ 

PERIOD 

(SEC) 

/ 

AMPLITUDE 

(DEG) 

1 

295-17:59:30 

295-18:02:31 

181 . 

9050 

323 

0.5 

/ 

5 

/ 

1.25 







1.0 

/ 

5 

/ 

2.5 







2.0 

/ 

5 

/ 

5.0 




: 



3.0 

/ 

5' 

/ 

7.5 

2 

295-18:10:15 

, 295-18:11:09 

53 • 

2650 

323 

Moderate Size Random Input 










3 

295-18:12:25 

295-18:13:21 

56 

2800 

323 

0.5 

/ 

10 

/ 

2.5 




• 


• 

1.5 

/ 

10 , 

/ 

■ 7.5 

4 

295-20:21:01 

295-20:22:05 

64 

3200 

324 

1.5 

/ 

5 

/ 

3.75 



- 


■ 


2.0 

/ 

5 

/ 

5.0 

5 

295:20:22:05 

295-20:23:04 

59 

2950 

324 

0.25 

/ 

50 

/ 

6.25 

6 

295:20:24:00 

295-20:25 

60 

3000 

324 

Step Inputs 



















Table 3.2 
Sensor List 


SENSOR 

ENGINE 

UNITS 

Throttle Coranand 

Left 

Deg 

Throttle Pickoff 

Left 

Counts 

Fuel Flow 

Left 

PPH 

Compressor Speed 

Left 

RPM 

Fan Speed 

Left 

RPM 

Exhaust Pressure 

Left 

Counts 

Duct Pressure 

Left 

Counts 

Throttle Command 

Right 

Deg 

Throttle Pickoff 

Right 

Counts 

Fuel Flow 

Right 

PPH 

Compressor Speed' 

Ri ght 

RPM 

Fan Speed 

Right 

RPM 

Exhaust Pressure' 

Right 

Counts 

Duct Pressure 

Right 

Counts 
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The flight data was recorded on an in-flight data 
acquisition unit and the digital tape processed to provide- 
scaling, units and conversion for use on the CDC 6600 
at NASA Ames Research Center. 

3.2 INITIAL PARAMETER ESTIMATES 

Classical techniques were used to calculate initial para- 
meter estimates using plots obtained from the data. Four areas 
were investigated using left engine data: (a]| steady state 
characteristics, Cb) hysteresis, (c) rate limiting and sym- 
metry, and (d) f an/corapressor dynamics. These are discussed 
in the subsections below. 

3.2.1" Steady State Operation 

An accurate model of Spey operation must match engine static 
behavior in this configuration. The engine model described 
in Chapter II provides a linear operating line for the state 
and control variables. Figure 3.1 is a compilation of measured 
engine response at steady power points. It will be noted that 
the approximation of linear spool speed to throttle response 
above about 5 degrees throttle is quite good. However, thrust 
and airflow are nonlinear functions of throttle. Since these 
are the primary model outputs, a tabular relationship can be 
developed for these variables as a function of rotor speed 
or throttle. Table 3,3 shows these functions tabulated against 
compressor speed. The cold and hot thrust figures refer to 
the thrust developed by the exhaust flow and duct flow of the 
engine at an assumed equilibrium rotor speed match. During 
transients., rotor speeds are not matched and an alternate 
formulation of these -.curves must be developed. 

An approximate set of curves for the thrust and airflow 
variables can be developed if it is assumed that fan airflow 
and core airflow are dependent only on the rotational speeds 
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Table 3.3 


Spey Engine Outputs as a Function of 
Steady State Corrected Engine Speed 


CORRECTED HOT STREAM THRUST T^/6 (Ibs^) 


t SPEED 

% SPEED 

INDICATED AIRSPEED - KNOTS ‘ 1 


N„//e 

V = 0 

V = 60 

V = 250 

0 

0 

0 

0 

0 

55.0 

84 

1023.2 

1092.1 

2102. Q 

70.5 

89,5 

■ 2402.7 

2456.7 

3230.0 

78.9 

92.5 

'371 4. 2 

3763.2 

4476.3 

86.0 

95 

4834.6 

4876.5 

5465.0 

94,4 

98 

■5867.0 

5902.6 

6411.0 

109.9 

103.515 

7013.5 

7042.1 

7446.7 

CORRECTEl 

) COLD STREAM THRUST - T^/6 (Ibs^) 

% SPEED 

% SPEED 

INDICATED AIRSPE 

ED - KNOTS 

Nj 


V = 0 

V = 60 

V = 250 

0 

0 

0 

• 0 

0 

55.0 

84 

796.7 

783.7 

914.6 

70.5 

89.5 

1660.2 

1666.7 

1807.2 

78.9 

92.5 

2494.5 

2489.6 

2668.4 

86.0 

95 

2977.9 

2954.4 

3110.4 

94.4 

98 

3301.0 

3289.3 

3442.9 

109.9 

103.515 

3532.3 

3554.1 

3682.4 





Table 3,3 CConcluded) 


% SPEED 

\/ye ' 

% SPEED 
N^/v/0 

CORRECTED ENGINE 
AIRFLOW 

m,/6/(Slugs/Sec.) 

a 

0 

0 

0 ' 

55.0 

84. 

2.514 

70.5 

89.5 

4,003 

78.9 

• 92,5 

4.964 

86.0 

95. 

5.507 

94.4 

-98. 

5.935 

109.9 

103.515 

6.488 
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of the compress ion elements. This is normally a more accurate 
assumption than the constant by-pass' ratio assumed in Table 
3.3, In this case, hot thrust is a function of the compressor 
speed only. Cold thrust and net airflow are determined by 
the fan speed. During transient motion each speed can vary 
independently . 

Data is required to determine the equilibrium rotor speed 
match. This data was measured from the flight tape when the 
engine had settled at a constant power. Figure 3.2 shows this 
match and verifies the assumed linearity of the match in this 
power range. Figure 3.2 and Table 3.3 can be used 
to determine a new functional representation of cold 
thrust and mass flow plotted against fan speed. Then, during 
simulation, nonequilibrium fan ..speeds are used to calculate 
instantaneous cold thrust and mass flow for the engines. 

Similar steady operating line data can be derived from 
the flight records for the intermediate simulation variables 
vis- fuel flow, throttle angle, and pressures as a function of 
the states of the simulation. Relationships derived from 
the data are shown in Figure 3.3 for fuel 'flow and throttle 
angle as a function of rotor speed. This line should be 
compared to Figure 2.14 which is the predicted linear relation- 
ship based on the governor characteristics. The operating 
assumptions of the linear operating linear hold quite well 
for this particular flight regime. Steady state linear 
operating equations are listed below: 

0 = 0.7212 05.J. -.31.5°)- (N^j - 100%) 

= 2.81 Njj - 181% 

='7.13 - 613% 

PTD = 1,32 - 32% 

PT6 = 4.0 Njj - 300% 

where the 100% values of the variables are listed in Table 3.4. 
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20.4 21.8 23.2 24.6 26.0 27.3 28.7 30.1 

(DEG) 


' 31.5 


Figure 3.2 Rotor Speed-Match (AWJSTOLRA 
Left Engine) Flight 323/324 




* 1% = 50 pph 


Figure 3.3 Fuel Flow Steady State Data 
(AWJSTOLRA Left Engine) 
Flight 




Table 3.4 

100% Values o£ Engine Variables 


VARIABLE 

— 

VALUE 

UNITS 

«T . 

31,5 

' Degrees 

% 

-957 

Counts 

“f 

5000 

pph 

"h 

12135 

RPM 

\ 

10075 

RPM 

PT6 

3775 

Counts 

PTD 

5036 

Counts 


It should be noted that the operating line linearity of 
engine variables, e.g.' the rotor speed match is predicted 
from previous engine experience hear full power. Governor 
linearity is predicted from the functional characteristics 
of the centrifugal counter-weight forces and the square root 
pressure drop - mass relationship used in the metering orifice. 
The latter relationship can be used to derive other character- 
istics of this control response as is discussed in Section 
3.2.5. 

3.2-. 2 Throttle Hysteresis 

Classical methods can be used to approximate the hyster- 
esis width since data containing periodic excitation is 
available. The ramp inputs are sufficiently close to sinu- 
soids that phase lag measurements can be used as a good 
measure of the hysteresis width. The sinusoidal describing 

function for hysteresis is shown in Figure- 3.4. Considering 
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Input 


Nonitnear Function 


Output 



0 0.1- 0.2 0.3 0.4 0.5 0.6 0.7 03 0.9 1.0 

atA 



Figure 3.4 Hysteresis Describing Function 
Characteristic 
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only the phase shift characteristics, measured values . 
can be related to the hysteresis width and input amplitude, 

The characteristic is independent of frequency which allows 
run-to-run comparison of results. The phase shift measurement 
allows a far more accurate estimate of the hysteresis width 
than amplitude roll-off when graphical data is plotted. An 
example of the phase shift measurement is shown in Figure 3.5. 
Hysteresis points were taken at several values of input magni- 
tude and frequencies. These points are shown in Figure 3.6 
with a least squares fit of the hysteresis characteristic: 

({) = fjjCa/A) + 4)^ 

where f^^ is shown in Figure 3.4, is a residual phase shift 

unexplained by hysteresis, and A is the amplitude of the periodic 
throttle inputs. 

The throttle curve in Figure 3.6 shows the phase- shifts , 
between the commanded throttle input and the throttle pickoff. 
Since no frequency dependent lags are present, all data points 
can be overplotted. The fit shows the hysteresis width to be 
about 0.48 degrees with only 2 ° of phase shift unexplained by 
the hysteresis. This residual phase is within the accuracy bounds 
of the data so it can be concluded that a simple hysteresis is an 
excellent model of the throttle linkage. 

Phase measurements were made between the high rotor speed 
response and both commanded and measured throttle inputs at 
a single frequency (o) = 1.26 rad/sec). The amplitude dependent 
portion of the response can be modeled as an additional hysteresis 
within the fuel control. The curve for the governor response 
indicates hysteresis width of about 0.19 degrees- This es'timate 
may be less accurate due to honlinearities in the governor 
response which will tend to cause amplitude dependent phase 
shifts. 

The hysteresis figures calculated above can be used to 
correct the phase response of the high rotor speed to throttle 
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(j)^ = Phase. lag between 6^ and 6.^. 


(|>u. = Phase lag between 6 and N,, 

n C ^ n 

(j)|^ = Phase lag between 6^ and 







THROTTLE AMPLITUDE 
(DEG) ■ 

A 


Figure 3.6 Identification of Hysteresis from 

Phase Lag 



inputs. This is shown in Figure 3-7 for various frequencies. 

A fit of a ‘first order lag phase response is shown which mini- 
mizes the sum squared error of the corrected points. The time 
lag associated with the fitted phase characteristic is 0.9 sec 
which should he close to the high rotor time constant (closed- 
loop) . 

3.2.3 Rate Limiting 


The governor operation described in Chapter II indicates 
that throttle inputs act as rate commands to the pressure drop 
governor. Since the inputs during the flight test .were throttle 
ramps, the rotor speed rates which are attained will reflect 
rate limiting by the governor if it exists.. Rates were ■ 
measured for various acceleration/deceleration magnitudes 
in the flight records. Rotor speed rates are plotted against 
throttle rates in Figure 3.8. The test run data was augmented 
by miscellaneous points which occurred when ramp inputs were 
generated as commands during the normal cockpit operation. 

For a linear operating line, and no rate limiting by the con- 
trol, the RPM response rate should be determined by the steady 
state operating line as is shown below: 

(Steady State) (N^^ - 100%) = 0.7212 (6^ - 31.5) 

CSteady Rate) = 0.7212 6.^, 

This line is drawn on the figure. The results show good cor- 
respondence to the predicted operation of the control. For 
throttle rates less than 5.5 deg/sec during accelerations or 
decelerations , the governor characteristic is linear. Rate 
limiting occurs at higher throttle input rates due to the 
saturation of the throttle sleeve on the stops. It appears 
that the acceleration and deceleration stops are not sym- 
metrically adjusted so that large engine declerations are 
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FREQUENCY (RAD/SEC) 


0.05 0.1 0.5 1.0 2.0 



Figure 3.7 Throttle to Rotor Speed Phase Difference (Data Corrected 

for Linkage Hysteresis) AWJSTOLRA Left Engine, Flight 323/324 
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Figure 3.8 Fuel Control Rate Characteristic 
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somewhat slower. For maneuvers during test sequences, the 
governor does not rate saturate and this effect can not be 
precisely identified. 


3.2.4 Fan and Compressor Dynamics 


In the development of the dynamics of the engine in 
Chapter XI, it was noted that -the engine rotors have two modes 
of response, namely, a unison mode and a fan rotor relaxation 
mode. The magnitude of these time constants can be approxi- 
mated from the flight data by examining the rotor phase relation- 
ship at various frequencies. This phase , angle is not dependent 
on the characteristics of the fuel control. It is assinted that 
fuel flow primarily enters the compressor speed equation, then 
the transfer function bet^^een the rotor speeds can be written 
as follows: 


Nh S*1 

where Tjj and are the time constants of the two* modes. 

The phase differences were measured at three forcing frequencies 
and are plotted in Figure 3.9. A function representing the 
first order phase difference was fit to the data, 

(j) - <j> CNjj) = tan ^ CtjjO)) - tan ^ . 

A time constant of 1.0 sec was determined for the compressor. 
This corresponds quite well with the fuel flow to rotor speed 
phase calculation corrected for hysteresis. The fan speed 
time- constant is predicted to be 2.6 sec. This is a rather 
large value for typical turbofan configurations without the duct 
bleed. Since only three points were used, the accuracy figures 
are unkno^rai and the estimate could be inaccurate . 
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FREQUENCY (RAD/SEC) 



Figtire 3.9 Spool Dynamics Identification (from Fan and Compressor Phase 
Difference) AWJSTOLRA Left Engine - Flight 323/324 



3,2.5 Control Parameters 


The control model developed for identification has a number 
of static elements which can be estimated from steady state and 
dynamic responses of the controlled engine. These calculations 
are discussed below. The assinnption is made that the null 
position in the speed governor is independent of engine speed. 
This is nearly true for moderate power changes . In this case,, 
the governor error, e, can be written from the steady state 
operating line as follows: 

N - 0.7212 5^ - 77.28 = -e C% RPM) 

The governor error is zero in steady state; thus, the pressure' 
drop governor must cancel the engine dynamics to form a pure 
integrator Cor rate servo) , This implies that the engine fuel 
to rotor speed gain, c^, is inversely balanced by the pressure 
drop regulator characteristic rotor speed to fuel flow Cpressure 
idrop) gain. This situation is represented in the model. The 
value of C£ is determined from the operating line relation- 
ship. 


= 7.13N - 613% 
or 

c^ = 7.13 (% F.F./% RPM) 

The remaining control parameters are the governor forward 
loop gain and nonlinearity, a. The control law for the engine 
and pressure drop regulator can be written as follows (tceglec- 
ting the nonlinearity) 




A A 


^F^AF 


CO. 72. 2 6^ - 


77.2,8) 
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The time constant of the closed- loop system is t = „ 

L»L . n r Ar 

and a reasonable guess of the value for K. is. 1,0 sec"l, 
The' nonlinearity , a, which results in slower accelerations 
than decelerations is difficult' to estimate. A nominal value 
of zero is chosen for the initial estimation runs. 

3.2.6 Summary of Initial Parameter Es'-tim'ates 


The discussion in Section 3.2 has provided the basis for 
initial parameter estimates for the engine/control model 
developed in Chapter II. These estimates are summarized in Table. 
3.5 and identified with the appropriate computer program symbology 
Csee Chapter Y for estimation results) . Initial estimates for 
scale factor and instrumentation biases are also included. 

These were measured directly from the flight data. 
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Table 3.5 

Initial Parameter Estimates 


PARAMETER 

SYMBOL 

INITIAL 

ESTIMATE 

UNITS 

Throttle Hysteresis 

DLINKl 

0.48 

Deg 

Control Hysteresis 

DLINK2 

0.19 

Deg 

Throttle Pickoff Scale Factor 

BPOTl 

-25 

Deg/Counts 

Throttle Pickoff Bias 

BPOT0 

950 

Counts 

Throttle/RPM Gain 

KTH 

0.7212 

Deg/% RPM ‘ 

Area to Fuel Gain 

KAF 

1.0 

Sec“^ 

Governor Nonlinearity 

ALPHA 

0.0 

(% RPM)"^ 

Flowmeter Damping Ratio 

ZETA 

0.2 

— 

Flowmeter Natural Frequency 

OMN 

13 

Sec"^ 

Compressor Lag 

CLAG 

1.0 

Sec 

Fan Lag • ' 

FLAG 

2.6 

Sec. : 

Fuel/RPM Gain 

CF 

7.13 

% RPM/% PPH 

Fuel to Fan Coupling 

CFL 

0.0 

% RPM/% PPH 

Fan Speed Bias 

XNLBIS 

0 

% RPM 

Tailpipe Pressure Lag 

TLAG 

0.0 

Sec 

Duct Pressure Lag 

DLAG 

0.0 

Sec 

Tailpipe Pressure Scale Factor 

TGAIN 

4.0 

%/% RPM 

Duct Pressure Scale Factor 

D6AIN 

1.32 

%/% RPM 

Tailpipe Pressure Bias 

PT6BIS 

-300 ■ 

% 

Duct Pressure Bias 

PTDBIS 

-32 

% 
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CHAPTER IV 

PARAMETER IDENTIFICATION PROGRAM 


A program was developed to process the flight acquired 
data and estimate parameters of the model. The details of 
the software are presented in this chapter. A program listing 
is included in Appendix A. 

4.1 THEORETICAL BACKGROUND 

The Spey engine model developed in Chapter II is used 
as the foundation of the parameter estimation procedure. This 
model is sho™ in Figure 4.1. The model' equations can be re- 
presented as follows: 

X = f(x, u, 0) (4.1^ 

= hCxj^, Uj,, 9) + C^.2) 

where the x is the engine state vector, y is the measure- 
ment vector, f and h(.,.,.) are the engine/control 
dynamic equations and measurement distribution matrix,- respec- 
tively. Table 4.1 lists the element of the state, output, 
and control vectors. The measurement noise is assumed to have 
a Gaussian distribution with known covariance. 

The maximum likelinood parameter estimates are obtained 
by minimizing the log-likelihood function for the parameter 
estimate given the data. A complete discussion of the theory 
is given in Ref. 6. 

J = min z (yCi)-(yj)V^C7Ci)-y(i)} C'l-S) 

'' i=l 

e 

2 2 

where R = diagCa^ >•••} ^ covariance matrix of the 

independent measurement errors, and' 

y(i) = hci(i), uCD, 9) (4.4) 
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Table 4.1 

Engine State, Output and Control Vectors 



SYMBOL 

STATE 


High Rotor Speed 

XNH 

Low Rotor Speed 

XNL 

Tailpipe Pressure 

PT6 

Duct Pressure 

PTD 

CONTROL VECTOR 


Throttle Position 

DELTT 

MEASUREMENT VECTOR 


Throttle Pickoff 

DELTPO 

Fuel Flow 

WFMSN 

High Rotor Speed 

XNHMSN . 

Low Rotor Speed 

XNLMSN 

Tailpipe Pressure 

PT6MSN 

Duct Pressure 

PTDMSN 


subject to the differential -constraint 

X = f (x, u, 0) . (4.5) 

For the Spey engine, Eq. (4.5j is discretized internally 
in the simulation program to form the discrete constraint equa- 
tions : 

x(i+l) = <j)x(i) + ruCi). (4.6) 

The minimization in Eq. (4.3) is carried out using either a Quasi- 
Newton method or Levenberg-Marguart. procedure. The. Levenberg 
Marquart algorithm provides a stable search for the minima near 
a priori estimates. The narameter step for this algorithm is 
as follows; 

2-1 

6(i+i) - ecu -C^ XI) || (4.7) 

8 6 
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where the expression, 



is approximated from gradient calculations and the Marquart para- 
meter, X, is chosen based on minimization performance during 
the search. More detail on these algorithms is contained in 
Ref. 6. 

The parameter estimation procedure described produces the 
maximum likelihood estimates for a given data record assuming 
the correct model structure and no process disturbances. These 
assumptions are well satisfied in the engines. The model 
structure has been carefully studied and initially evaluated 
to assure a good a priori structure. 

An additional identification option is included in the 
program to facilitate the initial startup of the procedures. 
Convergence of the maximum likelihood procedure is slowed by 
poor initial estimates and by improperly excited modes. To 
improve initial parameter estimates before the full algorithm 
is run, a partial identification procedure can be used. The- 
program logic allows utilization of measured variables as 
inputs to model blocks in the simulation. In this way, smaller 
identification problems can be run to obtain accurate a priori 
guesses for the full system.' As an example, measured fan speed 
can be used as one input to the duct dynamics. Gain and time 
constants can be found which then are used as estmates for the 
full identification. For the full run, the derived fan speed 
is used as the input and the overall response is estimated. 

4.2 STATISTICAL PROPERTIES 

The identification algorithm generates an estimation of 
the statistics of the parameters. These figures may be used to 
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evaluate the parameter accuracies and combine estimates from 
different data records to calculate the best overall para- 
meter results. Matrix outputs from the program are discussed 
below. 

The cost function, J(0) is calculated from normalized 
measurement error covariances as follows : 

P R N , 

J = Z ^ 2 (4.8a) 

i=-l. ^i j = l 

^i^^^ = 7iCj) - C4,-8b) 


where R = Ot 
o 1 

Thus, the measurement accuracies are weighted with respect to 
the accuracy of the first instrument. The scaled information 
matrix is derived from the second derivative of J with respect 
to the parameters. This matrix is developed during the iter- 
ative solution of the minimization problem.. The scaled inform- 
ation matrix is defined as follows: 


H 



(4.9) 


The eigensystem (eigenvalues and eigenvectors) of H are 
calculated. The overall conditioning of H and linear combin- 
ations of parameters corresponding to small eigenvalues of H 
are diagnostic tools in evaluating the convergence. An esti- 
mate of the parameter accuracies can be determined from the 
dispersion matrix as follows: 

D = N* 

and C4.10) 

cov(00'^) • > D 


where the matrix D forms the Cramer-Rao bound on the esti-- 
mates if N* = 1. 
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The value of N* is chosen to weight the independence 
of the data. A good approximation to this value is the ratio 
of the number of points between uncorrelated residuals, or, 
if the residual correlation t'ime constant is approximately 


determined to be t 


COR" 


then N* can be chosen as follows: 


N* 


'^COR 

TNtqt 


(4.11) 


where T is the sampling interval and is the total 

number of data points. Values of N* are typically chosen 
between 1 and 20. 

The fractional variance is calculated from the unit dis- 
persion matrix (D for N* = 1) . This equation is 


= A* °ii (4.12) 

where ^0^ is the fractional confidence interval on the esti- 
mate of the parameter, 9^. ’ • 

The F- ratio of the parameter estimate indicates the signi 
ficance of the parameter to the model. The F-ratio is the 
inverse of the fractional variance, or. 


F = ( 4.135 

1 1 

Finally, the fractional gradient of the likelihood func- 
tion indicates the ’’closeness’’ of convergence of the algorithm 
This quantity is defined as follows : 


Gfl = 3J/J/ae/9 ■ ■ (4.14) 

The quantities listed above are output for each parameter 
estimated. An example of this output is given in Figure 4.2. 
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Figure 4.2 Example of Spey Program Output Statistics 



FlUAl, RkSULIS Of UPllMllAUON 

0UA81**NEWT0H method 


0\ 

4^ 


FARM NO, 

variable 

location 

value 

1 

2 

ALO 

it 

.ZSITT^i 

ALT 

12 

,tt3U42E-02 

3 

OPOTO 

14 

927,896 

4 

epOTl 

IS 

•*2«,6Ag2 

Final root 

MEAN square 

RESIOUALS 


ERROR NO, 

RESIDUAL 

locaiioh 

VALUE 

J 

ImreRR 

1. 

4,0U3b’ 

average RH3 ERROR 

U,0163S 


VielCHT 

UOOOOO 


final 8TATISTIC3 


residual 3UH OF SQUARES = 52«25,‘> 

NORH Of THE CHADltNT a ,812278 

NO, OF pUNCr evaluations « 112,000 

ESI NO, OF SION, OieiTS “ 7,fa30bT 

final pARQUARbT parameter « 0, 

number Of ITERATIONS ■ » 0, 


NORMAL THRMINATIDN 


Figure 4.2 (Concluded) 


ORIGEfAU PAGE IS 
OF VOOR QUALmr 



Typically, for identification runs on different data 
records, the parameter estimates will be different-. This' 
is primarily due to the stochastic nature of the inputs. 

In some situations, accurate values of the parameters may 
be obtained when an optimum input time history excites the 
system response in a favorable way. Other data records may 
contain information about the parameters but not of the same 
quality as optimum inputs. The results of the parameter 
estimation runs can be combined to form the optimum estimate 
for the overall series of data records by weighting the average 
by the predicted parameter covariance. 

The exact formula is given in terms of the information 
matrices, M^, for 'each of q runs: • 


0 



q 

l M. 



C4.15) 


1 - 1 i 

where M. .= rr* R H. and 0. is the parameter estimate for 

1 N? O ; 1 X 

the ith record. An approximate formula can be used with reason- 
able accuracy as follows: 


8 - 

1 



ir 





C4.16) 


where is the jth diagonal element of the i' dispersion 

j 3 

matrix. The above formula can be used to obtain weighted 
averages of the parameter values accurately reflecting the 
overall responses. 
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CHAPTER V 

RESULTS OF PARAMETER IDENTIFICATION 


The model developed in Chapter II was used as the founda- 
tion of a parameter estimation program. In this way, the code 
used for identification could be directly transferred to the 
AWJSTOLRA simulation without reformulation of the modeling 
equations. The parameter values would be the same as identi- 
fied from the data. 

The identification procedure is described in Chapter IV 
and a listing of the simulation/identification program is 
included in Appendix A. Preliminary parameter estimates were 
graphically calculated from the flight data as described in 
Chapter III, 

The simulation can be used to identify parameters in 
portions of the engine using measured or derived variables as 
input. The parameters which. are identified can be optionally 
selected. Table -5.1 shows the parameters of the engine simu- 
lation which are available to estimate. Variable name defini- 
tions associated i^ith block diagram variables (Figure 5.1) 
are shoivn in Table 5.2. The program was sequentially exer- 
cised to identify throttle parameters, engine pressure dynam- 
ics, engine dynamics, and engine/control dynamics. The 
identification setup used and the results are presented in 
the following subsections. 

5.1 THROTTLE IDENTIFICATION 

The throttle hysteresis value and pickoff bias and scale 
factor were identified. In addition, a power level dependent 
hysteresis width was hypothesized and added to the parameter 
list. Figure 5.2 shows the program set-up. Measured throttle 
commands are used as input, the throttle pick-off is commanded 
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Table 5.1 


SPEY engine simulation PARAMETERS 
AWJSTOL RESEARCH AIRCRAFT 
. SYSTEMS CONTROL? INC tVT) 

1801 PAGE mill ROAO 
PAUO ALTO? CALIFORNIA ■ 



* * 

* RIGHT .engine T 

* FLIGHT NUMBER 32«* 

* RUN NUMBER 3 A 

* STI 58000 * 

♦ ..*.* ♦ ♦ 

♦ throttle T control ♦ ENGINE ♦ DUCT/EXhAUST * 

♦ ♦ * * ■ ■ ■ * 

♦ T*YT»ie*******T + T*»*»*»*T!Ti**+*»^?:»»:T*?!*»Tt**+**»1‘**»»4**»»f »*?♦*♦»< 

♦OLINKja ,329645 ‘ *KTH » ,782<»«6 *CLAG s ,210000 ♦TLAG = ,SOOOOOE"02* 

♦0LINK2» ,lOOOOOE-Ol*NAF s 5,04400 *FLAG * ,403000 ♦TGAIN a 4,72200 ♦ 

♦ ALO = ,253775 *ALPHA * 0, 4CF * ,110000 • ♦PT6Bl3a-»l ,94100 * 

♦ALl 3 .437142E-02TCFCNTR3 ,110000 *AFL » .710000 *OLAG ,50(JOO0E-O2* 

♦AL2 3 0, aERSBISa ,435724 *CHL » 3,07300 *0GAIN = 1,13300 ♦ 

♦BPOTO 3 927,898 *KF3F b 1,00000 *XNLHI5» ,324000 ♦PT0BIS343, 93000 * 

♦ BPOTI »*24,569H TKFaiA3* 20,9200 .♦CLAGO » ,909156 *TLAGO 3 ,1831-56E<.0l* 

♦ZETA s ,990000 TOM a 50,0000 *CLAGG a ,908436E»01*TUAGG = ,901664 * 

♦OH. 3 50,0000 ♦SIG a-,202oaoE?>01*FLAGD » ,951564' ♦DLAGO a ,1831362-019^ 
♦KAt 3 0, *KA2 a 0, .*FLAGG » ,4a4l65E-0l90LAGG a ,961684 ♦ 


♦♦♦♦♦»?♦♦♦♦♦ ♦♦♦♦♦♦ ♦♦♦♦♦:»♦*♦♦♦♦♦♦♦♦■?♦*♦♦♦♦ 

♦ FULL POWER POINT * 

♦ ■ ' ■■ ■ ■ * 

♦ throttle POSITION a 31,30 DEGREES * 

♦ throttle PICKQFF a 976,0 COUNTS' * 

♦ MAIN fuel flow 3 5750, LBM/H0UR9 

♦ COMPRESSOR speed = ,1335E*03 REV/HIN ♦ 

♦ FAN 3FEE0 a ,i008£+05 REV/HIN ♦ 

f AUGHENTOR pressures 50,00 LBF/3QIN» 

♦ OUCT pressure s 50,00 LBF/SQIN* 

♦ ♦ 


♦♦♦♦♦♦♦♦a ♦♦♦♦♦♦^♦♦9 *♦♦♦♦=?*» ♦♦♦♦»»♦' 

♦ optimization parameters ♦ 

♦ ♦ 

* NO, OF measurements a 6 ♦ 

* NO, OF parameters 'a 4 ♦ 

* NO, OF ERROR TERM-S a ' 1 ~ ' ♦ 

* NO, OF DATA POINTS a 3250 » 

* NO, OF FIRST POINT * 1 • * 

* NO, OF increment = 1 ♦ 

* DATA PRINT INTEHVAL9 0 ♦ 

* ■ _ ♦ 

♦ CONSTANT print FLAGa ' q « 

♦ epsilon" a ,l0000£i*039 

♦ DELTA a .lOOOOE-03* 

♦ MAX FUNCTION tVALS a 300 ♦ 

♦ NO.OF'SICNIF OIGItSa ‘ 4 ♦ 

♦»♦♦♦»♦♦*»♦♦♦♦**♦♦ ♦♦♦»♦»♦♦? ,!»?♦*»♦ 
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Figure 5.1 Proposed Spey Engine 



ion Model Showing Program Symbology 





































Table S.2 

Spey Parameter Index 


PARAMETER 

INDEX 

SYMBOL 

NAME 

11 

ALO 

Constant throttle hysteresis 

12 

ALl 

Linear throttle hysteresis 

13 

. AL2 

Quadratic throttle hysteresis 

14 

BPOT0 

Throttle pickoff bias 

15 

BPOTl 

Throttle pickoff scale factor 

16 

DLINK2 

Internal control hysteresis 

17 

KTH 

Throttle to RPM transfer factor 

18 

ALPHA 

Throttle nonlinearity constant 

19 

KAF 

Area to fuel gain 

20 

CFCNTR 

Pressure drop governor feedback 

21 

ZETA 

Fuel flowmeter damping ratio 

23 

OMN 

Fuel flowmeter natural frequency 

26 

CLAG 

Compressor speed time constant 

29 

CF 

Fuel to RPM gain 

30 

AFL 

Fuel to low rotor speed gain 

31 

CHL 

High rotor to low rotor gain 

32 

FLAG 

’ Fan speed time constant 

•35 

TLAG 

Tail pipe pressure time constant 

38 

TGAIN 

Tail pipe pressure to RPM gain 

39 

DLAG 

Duct pressure lag 

42 

D6AIN- 

Duct pressure to RPM gain 

43 

KAl 

Flowmeter estimator gain 

44 

KA2 

Flowmeter estimator gain 

46 

EPSBIS 

■Throttle bias 

47 

KFSF 

Flowmeter scale factor 

48 

KFBIAS 

Flowmeter bias 

50 

PTDBIS 

Duct pressure bias 

51 

PT6BIS 

Tail pipe pressure bias 

52 

XNLBIS 

Fan speed pressure bias 
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Figure 5.2 Setup for Throttle Hysteresis Identification 


with the derived pick-off signal and parameters adjusted to 
maximize the fit. Program convergence was quite consistent 
on a run-to-run basis. .The fit results for two runs are 
shown in Table 5.3. An examination of the fit errors indicate 
that the right engine modeling accuracy is significantly poorer 
than the left engine. The data shows a large amount of noise 
on both the throttle command and pick- off signal relative to 
the left engine. This noise is averaged by the identification 
procedure and does not significantly affect the accuracy of the 
resulting parameter estimates. Table 5.4 shows that the power 
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Table 5.3 

Throttle Input Fits 




NO. OF 
POINTS 

RMS THRERR 

RMS THRERR 

RUN 

RECORD 

LEFT 

(COUNTS) 

RIGHT 

(COUNTS) 

1 

1 

6000 

1.67 

6.10 

2 

5 

3250 

1.49 

3.02 


level dependent hysteresis width is not a significant parameter 
in the model. This may be inferred by the small F-ratio and 
large fractional standard deviation. Table 5.5 shows the 
average values of the final model for the throttle hysteresis 
and picic-off. 

5.2 IDENTIFICATION OF PRESSURE 'bYNAMICS 

The tailpipe and duct dynamics were modeled as pressure 
lags. It was expected that exhaust pressure lag would be 
insignificant. The duct lag was more uncertain. An attempt 
was made to identify the terms in a linear model of the flow 
as shown in Figure 5.3. The rotor speed measurements are used 
as input quantities. The error term is calculated from the 
measured and derived pressure. The terras modeled are the time 
lag, flow gain, and bias. Sensor scale factors are lumped 
with the flow gain. Input and output sensor biases are com- 
bined in the bias term. The flow lag and any sensor response 
lags are modeled in the time constant term. Typical transient 
pressure sensor response times are less than 10 msec and the 
rotor speed lag is less than this. Any identified lag 
should be associated with the flow dynamics in the exit 
volume. 
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Table 5.4 


Parameter Estimation Results for Throttle Hysteresis Using Flight Data 

of Records 1 and 5 


LEFT ENGINE 


; PARF 

1 NO, 

VARIABU£ 

VALUE 

fract,sto,dev. 

F RATIO 

FRACT.SENSn, 

1 

i 

Abo 

,251611 

,24965oE'»0 Ji 

1604,49 

,246065E«05 


2 

ALl 

I42ai9£«03 

■,2,20025 

,206565. 

t,146542Ei07 

© 

3 

bPOTO 

976,087 

. .172528E-03 

,335955Et08 

«,135067E-»03 

4 

BPOn 

*,26,2088 


,115ai9Et08 

•7S3474Ep04 


1 NO, 

varuhue 

VALUE • 

FRACT,SlD,ptV, 

f RATIO 

FHACf ,SENSn, 


\ 

ALO 

,253676 

, 1-42065 

48,9945 

«,'784314E«.0S 

(?) 

2 

AU 

, 199695E«03 

10,4109 

,92l2t0E-02 

w,l735e3E»06 


3 

•DPOTO 

986,654 

,173938E«03 

,33O532EfO0 

»,20060OE^O3 


4 

bpqt J 

^•26,3713 

■<(.,285379E«03 

,12ar88E+08 

,223072E»03 

...RIGHT ENGINE ' 

PARH 

1 NO. 

VARiAflbE 

VALUE 

FRACT,STP,DEV, 

F RATIO 

FRACT,3fcN3lT, 


i 

AlO 

,250647 

,163424 

37,4430 

,lll249E-05 

© 

2 

AU 

, 159342E*»02 

1,26673 

,623203 

,5475S7Ei>07 

3 

BPOTO 

936i661 

• .156515E-03 

,4082l6Et08 

w,415084E*p03 


4 

bpoti 

«24,9t0a 

'•,243666E<703 

\ 

,166l50Et06 

,267171Eo>03 

PARM 

NO, 

variable 

VALUE 

FRACl .STD.OEV, 

F RATIO 

FRACT.SENSIT, 


1 

ALO 

,253775 

,988602E«02 

10227,8 / 

og 184403E»06 

(?) 

2 

ALl 

,437H2E*'02 

,3329B9t.«0l 

901,064 / 

s.SyByOBfcoOT 


3 

BPOIO 

927,090 

•,1508l2E»O3 

,396492E*0b 

,226552E»04 


4 

BPOTI 

»24,5692 

«,252S30L«03 

,1568t0E+0B 

o,14515iE»04 








Table 5.5 ' 

Combined Throttle Estimates 



RUN 1 * 

EST/FSD (%) 

RUN 2- 

EST/FSD (%)* 

AVERAGE 

EST 

UNITS 

LEFT 





ALO 

0.252/2.5 

0.254/14 

0.252 

Deg 

BPOTO 

977/. 17 

987/. 17 

982 

Counts 

BPOTl 

-26. 2/. 27 

-26. 4/. 29 

. -26.3 

Counts/Deg 

RIGHT ■ 





ALO- 

0.251/16 

0.254/1.0 

0.254 

Deg 

BPOTO 

937/. 16 

928/. 16 

942 

Counts 

BPOTl 

"-24. 9/. 24 

-24. 6/. 25 

-24.75 

Counts/Deg 


* EST/FSD = Estimate/Fract, Std. Dev^ 


PT6BIS 
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An estimation run on record one was made. Results are 
presented in Table 5.6. The fits are shoivn in Appendix C. It 
is apparent that the linear model of the pressure is not completely 
satisfactory in measuring large pressure responses. This is prob- 
ably due to a combination of sensor and pressure response nonlin- 
earity. The estimated pressure gains and bias are best fits to 
the particular data record. The exhaust lag is identified as 
zero reenforcing the intuitive conclusion. The duct lag time 
constant is 0.049 sec. This figure represents an insignificant 
value relative to the engine dynamics. It is consistent with 
the lumped volume time constant of the flow, 

5.3 I DENT IP I CAT I ON OF ENGINE PARAMETERS. 

The parameters determining the dynamic* response of the 
engine to fuel flow inputs were identified. The program setup 
is shown in Figure 5.4, The measured fuel flow is the input used. 
Error signals are derived from rotor speed measurements and 
pressure measurements. Pressure gains and scale factors are 
introduced a*s free parameters 'in addition to the linear model 
engine equations. Data records representing triangle wave 
and doublet inputs were used. 

Table 5.7 lists the fit errors for identification of left 
and right engine parameters for the three sets of data. Rep- 
resentative time history matches are included in Appendix C. 

The results of the estimation runs are shown in Tables 5,8 and 
5.9 for the left and right engines, respectively . In general, 
the engine variables are accurately and consistently estimated. 
Intermediate values such as pressure gains and bias are not 
as consistent. These terms have been more accurately identi- 
fied in setup in Section 5.2. The inconsistencies arise be- 
cause of a nearly unidentifiable' combination of bias and 
scale factor on the instruments and pressure response nonlinearity. 
This does not affect the accuracy of the engine variables. 
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Table 5.6 


Identification of Pressure Dynamics 
(Record No. 1, 6000 Data Points) 


PARAMETER 

LEFT ENGINE 

RIGHT ENGINE 

UNITS 

PARAMETER 

FRACT. 
STD. DEV. 

PARAMETER 

FRACT. 
STD. DEV. 

TLAG 

0.0* 

- — 

0.0* 

— 

Sec 

TGAIN 

4.10 

0.008 

4.35 

0.03 

% PRESS/% RPM 

DLAG 

0.05 


0.05 


Sec 

DGAIN 

1-.36 

0.003 

1.40 

0.07 

% PRESS/% RPM 

PTDBIS 

-3.02 

0.01 

-3.53 

0.26 

% PRESS 

PT6BIS 

-5.41 

0.03 

-1.82 

0.40 

% PRESS - 


* Not- identifiable. 

** Estimate- not reliable because of large sampling ^period relative to 
lag. 
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Figure 5.4 Engine Parameter Identification 









Table 5.7 

Fit Errors for Engine Parameter 
Identification Runs 


LEFT ENGINE 





AVERAGE 

DATA 




ERROR% 

RECORD 

ERROR NO, 

residual 

location value 


1 

XHSERR 

3 

,148043 

1 

2 

XUSERR 

4 

,192332 


3 

PT6ERR 

. 5 

,378965 


4 

ptderr 

6 

,456491 


error 'no. 

RESIDUAL 

location value 


1 

XHSERR 

• 3 

,260625 

3 

2 

xlserr 

4 

,276202 


3 

PT6ERR 

5 

.281031 


4 

ptderr 

6 

,331822 


ERROR NO, 

RESIDUAL 

LOCAriON VALUE 


i 

XHSERR 

3 

,369610 

7 

2 

XUSERR 

■ 4 

• ,565856 


3 

PT6EHK 

S 

,417991 


4 

ptderr 

6 

,79l'49i 

RIGHT ENGINE 





AVERAGE 

DATA 




ERROR % 

RECORD 

ERROR NO, 

residual location 

VALUE 


1 

XHSERR 

3 

,121402 

1 

2 

xlserr 

4 

,991445E«»01 


3 

PT6FRR 

5 

,258342 


4 

PTDERR 

6 

,282912 


error no. 

residual location 

VALUE 


i 

xhsehr 

3 

1,74624 

3 

2 

xlserr 

4 

1,73914 


3 

PTftERR 

5 

,291593 


4 

PTDERR 

b 

,311548 


ERROR NO, 

residual location 

VALUE 


i 

xhseRR 

i 

,265544 

7 

2 

xlserr 

4 

,364427 


3 

PT6FRR 

5 

,279910 


4 

ptderr 

b 

,467315 
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Table 5.8 


Engine 


Parameter Identification Results 
for the Left Engine 




DATA RECORD 

1 

' 

3000 POINTS 


60 SECONDS 

paRm no. 

VARIA3LE 

VALUE 

FRACT.STO.DEV, 

F RATIO 

fract.sensit. 

t 

CLAS 

♦,24552a 

,183116E-01 

2982,27 

-,2669025*02 

2 

CF 

,112700 

,22372SE«02 

199789, 

-.963444 

3 

AFl 

,608423 

,6517«1E-01 

235,423 

-,8053315*02 

4 

CHI, 

2,70506 

,ll2671t-01 

7877,25 

-,927046 

S 

flag 

,395805 

.828014b-02 

14585,6 

,344291 

b 

TGAIN 

a,9Jl53 

,360138£*02 

77101,5 

,273931 

7 

OCAIN 

1,47099 

,l42866E*0i 

4899,40 

-*,515400 

8 

KFSIAS 

23.6544 

,394499£*02 

64255,2 

-,538333 

9 

PT08IS 

•3,19989 

•,540682E-01 

342,070 

-,126345 

10 

PT8HIS 

•1,75263 

",533703fc*0l 

351,075 

•,2451955*02 

U 

XNLB15 

,273644 

,214437E«01 

2174,70 

,427091 

1 - — - . 1 

DATA RECORD 3 


3000 POINTS 


60 SECONDS 

park no. 

VARIABLE 

- VALUE 

FRACT,STD,D£V, 

F RATIO 

FRACT.SENSIT,' 

1 

CLAG 

,209854 . 

.127346E-01 

6166,33 

.104873 

;2 

CF 

♦ 110279 

,345332E'-02 

83854,3 

,543096 

i 

AFL 

,709761 

,8345935*01 

143,566 

,7638385-Ql 

4 

CHL 

3,07253 

,6681i2E*02 

22402,8 

-.750232 

5 

flag 

,403171 

,3482505-01 

824,553 

.359266 

b 

TGAIN 

4,72230 

,3264365-02 

93843,4 

2,04080 

7 

dgain 

1.13533 

,644lSaEf02 

24099,9 

*,152601 

a 

KFBIAS 

20.9181 

,722847£.-0a 

19138,5 

.260554 

9 

PTOBIS 

•3,93060 

-,1605485-01 

3879,62 

-,216992 

to 

PT6BI3 

*1,94100 ■ 

•,376977£*01 

703,671 

,113070 

'll 

XNLBIS 

,324064. 

. ,23lll35»01 

1872,20 

,550566 

r 

DATA RECORD 7 


2000 POINTS 


40 SECONDS 

PARM NO. 

VARIABLE 

VALUE 

FRACT.STD.OEV, 

F RATIO 

fbact .sENsn. 

1 

CLAG 

,280539 

.174310E-01 

3291,23 

-,7iiaao5*0i 

2 

CF 

,132310 

,507232£,.02 

38867,5 

-19,0540 

3 

afl 

,648202 

,4744505*01 

444^241 

,33784lE*0l 

4 

CHL- 

3,55974 

,3012225-02 

110212. 

*6 ,65612 

5 

FLAG 

,471682 

,5506595-02 

32978,8 

2.23619 

6 

TGAIN 

4,03730 

.6373395-02 

24618,4 

-10,9313 

7 

DGAIN 

1,33361 

,7818485-02 

16358,9 

-2,90139 

8* 

XF81AS 

11,1947 

,3054455-01 

1071,33 

*3,43942 

9 

PTD81S 

•4,58556 

-.4273525-01 

.547,555 

*,466905 

to 

Ptaais 

*6,53689 

*,2376555-01 

1770,55 

•3,09991 

u 

XNL83S 

,387423 

,3299645-02 

91847,5 

2,25761 
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Table 5.9 


Engine Parameter Identification Results 
for the Right Engine 


DATA RECORD 1 



3000 POINTS 


60 SECONDS 

PASH NO, 

VARIABLE 

VALUE 

fract.std.oey. 

F ratio 

FRACT.SENSn. 

1 

CLAG 

,309677 

,269890E"01 

1372.86 

,497044E"01 

S 

CF 

,112470 

,3Sia39E-02 

81057,9 

3,7208/ 

3 

AFL 

,409645 

,500694E-01 

398,892 

-,390389E"0l 

4 

CHI. 

2,61253 

,880933E"02 

12885,9 

5,46837 

5 

FLAG 

,363951 

,115059E-01 

7553,68 

-1,96815 

b 

tgain 

5.17007' 

,59986aE"02 

62541.2 

"1.95474 

7 

OGAIN 

1,64233 

,156479£"01 

4084,03 

2,17827 

8 

KFBIAS 

24,6401 

,659776E"02 

22972,4 

1.61175 

9 

PTB81S 

•2.59336 

•.846524E-01 

139,547 

,366318 

to 

PT68IS 

1,96712 

,S622S3E-0l 

316.327 

,141740 

11 

XNLBlB 

.302576 

.105661E-01 

8957,20 

•1 ,96169 

L ■ - ... 1 

DATA RECORD 3 


3000 POINTS 


60 SECONDS 

PAJ?» »Ot 

variable 

VALUE 

FRACT.STD.DEV, 

F RATIO 

FRACT.SENSiT, 

1 

CLAG 

,319747 

,24943SE-01 

1606,61 

".123973E-03 

2 

CF 

,119920 

,338874E-02 

8708-1.0 

»,694889E"01 

3 

afl 

,7012/a 

,677632E-01 

217,777 

,190973Ei03 

U 

CHL 

3,15/78 

.441653E-02 

51266,9 

-,111738E-0l 

5 

flag 

,478460 

,680962£-02 

21565,3 

.120563E-01 

b 

TGAIN 

4,81038 

,344225E-02 

84394.9 

*,321819E"0I 

7 

OGAIN 

1,30134 

,938I84E-02 

11361.2 

•",a42486£i02 

8 

KF8IAS 

18,7169 

,867245E-02 

13295,9 

-,177729E"01 

9 

PTOSIS 

-10,0272 

-,124537E-01 

6447,71 

".329149E-02 

to 

ptaais 

"2,72262 

",a00297E»01 

624,074 

",413596E"0i 

11 

XNL31S 

,375867 

,261079E-02 

146709, 

,174602£"0t 

1 ■ ' " 1 

DATA RECORD 7 


2000 POINTS' 


40 SECONDS 

paRh No. 

variable 

VALUE 

FRACT,3TO.O£V, 

F RATIO 

FRACT.SEN3IT 

1 

CLAG 

,308398 

,l9I228£-0t 

2734,62 

,ll9934E-0 

2 

CF 

,155032 

,884564£"02 

12780,3 

•.633109E-0 

3 

AFL' 

", 161548 

",273005 

13,4171 

,4S5659E"0 

4 

CHL 

2,91548 

.963915E-02 

10762,7 

,202993 ■ 

5 

FLAG 

,249392 

,390653E-01 

655,267 

",S96218E"0 

b 

TGAIN 

3,40194 

,968594E-02 

10659,0 

",487642 

7 

OGAIN 

1,12496 

♦977269E-02 

1047.0,6 

.998672E-0 

a 

kfbias 

3.60752 

,130659 

53,5766 

«,523410£"0 

9 

PTOSIS 

"6,04426 

",3314S5E"01 

910,227 

,2768^0E"0 

to 

PT60I3 

"10,7760 

•»241091E"01 . 

1720,43 

".226384 

It 

xnlbis 

,379122 

,58a912£«02 

23833,7 

,406S96£"0- 
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Table S.IO shows the results averaged by the estimated 
variances. Accurate -estimates o£ rotor time constants and 
coupling were obtained from the data processing. The engine- 
to-engine variation in parameters is less than 10% except for 
the compressor time constant. The data indicates that the 
right engine time constant is 35% larger than the left, engine. 
Compressor and turbine build differences, seal leakage, bleed 
differences, deterioration, etc., are possible reasons for 
this type of variation. It should be observed that the fan 
and compressor response roots differ by about 50%. Modeling 
engine transients with a single time constant can lead to 
erroneous results as described in Chapter II. 

5.4 IDENTIFICATION OF CONTROL PARAMETERS' • 

There are three parameters in the Spey which determine 
the characteristics of the hydromechanical speed governor as 
modeled in the simulation. The setup for identification is 
shown in Figure 5.5^ The assiimptions made in the analysis of 
the governor produce this simple form. The results of the 
parameter estimation will verify these assumptions. 

The assumption that the positive feedback from the pres- 
sure drop governor just cancels the natural engine damping is 
an important simplification in the simulation. The rationale 
is briefly reviewed below. The transfer function from fuel to 
speed is a 'first order lag. If it is assumed that the pressure 
drop governor provides positive feedback with gain Cp, then 
the metering area to RPM transfer function is as follows: 

. Cp 

-!r “af 1 - Cp/c' 

where 

“^PD 1 - Cp/Cp 


TpjjS + 1 
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Table 5.10 

Parameter Estimates for Engines 


PARAMETER 

LEFT ENGINE 


r 

RIGHT ENGINE 

6 

r 


•RUN 1 

RUN 2 

RUN 3 

RUN 

1 

RUN 2 

RUN 3 

6 

r 

6 

r 

e 

r 

0 

, r 

0 

r 

6 

r 

CLAG 

.246 

.018 

.210 

.013 

.281 

.017 

.232 

iOI6 

.310 

.027 

.320 

.025 

.308 

.019 

.312 

.024 

Sec 

CF 

.113 

,002 

.110 

.003 

.133 

,005 

.114 

.003 

.112 

.004 

,119 

.003 

.155 

.009 

.118 

.005 

% RPH/5S FF 

AfL 

.606 

.065 

.710 

.083 

.648 

.047 

_.64i 

.065 

.410 

.05 

.701 

.068 

(-.16) 

C-27)* 

0.56 

.06 


CHL 

2.71 

.011 

3.07 

,007 

3.56 

,003 

3.39 

.007 

2.61 

.009 

3.16 

.004 

2.92 

.010 

3.02 

.008 

— 

FLAG 

.396 

,008 

.403 

.035 

.472 

,006 

.437 

.016 

.364 

.012 

.478 

.007 

.249 

.039 

.423 

.019 

Sec 

XNLBIS 

.274 

.021 

.324 

.023 

.388 

,003 

.382 

.016 

.303 

.011 

.376 

.003 

.379 

.006, 

.370 

.007 

% RPH 


* 

Not used In average 

0 a Parameter estimates 
r a Functional standard deviation 













Clearly, when Cp = Cp, the pressure drop governor has caused 
the engine to behave like a rate servo. The speed governor 
has a control as follows (assiime a linear metering character- 
istic) : 






and the closed- loop transfer function can be written as follows: 


where 


5N. 


H 


KthP 


^AF^F 


s + 1 


U = 


1 + 


^1 - Cp/Cg 
^AF^F 



which indicates that the parameter Cp i-s difficult to identify since 
Kpjj -proportionally determines the d.c. gain and Kaf proportionally 
determines the closed- loop time constant. The "identif iability” 
of Cp decreases as y approaches 1. There is a good justifi- 
tion that y is designed to be nearly 1. This can be seen 
from the metering area d.c. gain as follows: 


1 - Cp/Cp 

1 - Cp/C- + '^T 


(d.c.) 


It is desirable to hold the 'metering area at null since 
the' acceleration limits are determined by physically limiting 6A^. 
Null d.c. 5A^ is held if Cp =-Cp. This, assumption is made in 
the simulation. Mismatches in those two values cause the esti- 
mated value of K^p to be slightly biased, but this effect 
should be second order since y 1. 
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Using this assumption, K^P 


is identified with the. closed- 


loop response determination and with the operating (static) 

behavior. The values of are accurately determined by 

steady ground trim runs and verification -at selected in-flight 
points. This value is not estimated. 


The results of the estimation for the parameter, 
are shown in Table 5.11 for two data records. Appendix C 
contains sample plots of the time history match to the pre- 
dicted fuel flow. The loop time constant, 

calculated in the table. These values are reasonably close 
to those obtained from the graphical analysis. 


5 . 5 SUMMARY 

Component blocks of the Spey engine simulation have been 
used to estimate parameters from the flight data. Accuracy 
values have been calculated for each configuration. Several 
runs have been lised to generate each set of results. In 
general,, the run-to-run consistency has been excellent. Table ■ 
5.12 lists the final parameter estimates for each engine and 
the uncertainty limits associated with these estimates.. The 
confidence interval suggested for the parameters is calculated 
using the fractional standard deviation and a subjective 
estimate of the overall fit of the data estimation. Worst 
case values may lie with an interval two or three times larger 
than the suggested bounds. 
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Table 5.11a 

Identification Results for Engine Control, 
KAF and Parameters 


— 

RUN 

RECORD 

NO. OF 
POINTS 

ENGINE 

KAF 

FRACT. 
STD. DEV. 

1 

1 

3000 

Left 

■ 3.44 

,007 

2 

1 

6000 

Left 

2.76 

.008 

3 

7 

2000 

Left 

3.08 

. .012 

4 

3 

3000 

Right 

1.60 

.007 

5 

7 

2000 

R1 ght 

3.07 

.002 


Table 5.11b 
Averaged Values 



KAF 



{% AREA/% RPM) 


LEFT 

3.11 

0.65 

RIGHT 

2.96 

0.89 
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Table 5.12 

Final Spey Engine Parameter Estimates 


PARAMETER 

INDEX 

LEFT ENGINE 

RIGHT ENGINE 

ESTIJiATE 

UNCERTAINTY 

ESTIMATE 

UNCERTAINTY 

AL0 

U 

0.252 

+0,062 

0.254 

+0.060 

ALl 

12 

a/ 

1/ 

a/ 

a/ 

BPOT0 

14 

982 

+5.0 

942 

+5.0 

BPOTl 

15 

-26.3 

+0.22 

-24.8 

* +0.18 

DLINK2 

16 

a/ 

a/ 

y 

0.7212^ 

y 

KTB 

17 

0.7212^ 


ALPHA 

18 

Qe/ 

— 

0^ 


KAF 

19 

3.11 

+0.07 

2.96 

0.07 

CF 

20 

0 .II 4 S/ 

+0.001 

0.118^/ 

+0.002 

ZETA 

21 

0.2-'^ 

+0.2-'' 

0,2 

+0.2 

OHN 

23 

13^/ 


13^/ 

+3^/ 

CLAG 

26 

0.232 

+0.011 

0.312 

+0.022 

CFCMTR 

29 

0,114 

+0.001 

0.118 

+0.002 

AFL 

30 

0.64 

+0.13 

0.56 

+0.10 

CHL 

31 

3.39 

+0.07 

' 3.02 . 

+0.07 

FLAG 

32 

0.437 

+0,02 

0.423 

+0.02 

TLAG 

35 

0 

+0.001 

0 

+0.001 

■TGAIN 

38 . 

4.10 

j^.lO” 

4.35 

+0.40 

OLAG 

39 

0.05 

— 

0.05 


DGAIN 

42 

1.36 

+0.01 > 

1.40 

+0.3 

KFSF 

47 



1.0^'' 


KFBIAS 

48 

0.0^^ 

— 

0.0^'' 

— 

PTDBIS 

50 

-3.02 

+0.09 

-3.53 

+2.75 

PT6BIS 

51- 

-5.41 

+0.5 

-1.82 

+2,18 

XNL3IS 

52 

0.382 

+0.03 

0.370 

+0.03 

(sec) 


o'. 65 

+0.05 • 

0.S9 

+0.05 


— Not significant 

h/ 

— Estimate from inspection of data 
■^Defined as equal *to CF - See Chapter II 
-■^Derived from steady-state operating line 

fi/ . * 

— No significant estimate derived from data 
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CHAPTER VI 

SUMMARY AND CONCLUSIONS 


6 . 1 SUMMARY 

The Spey engines on the AWJSTOL research aircraft re- 
present a complex system of coupled dynamics. Flight data 
was acquired during controlled throttle tests in the powered 
lift approach configuration of the aircraft. A dynamical model 
of the' throttle linkage, control system, engines, and duct was 
developed from an analysis of the physical configuration of 
the engine. The flight data was utilized to develop initial 
estimates of parameter values using classical frequency re- 
sponse and describing function techniques. A system identi- 
fication program was usbd to process the flight and calculate 
estimates of the parameter value in the engine model. The 
resulting engine simulation should provide an accurate com- 
ponent repre-sentation for inclusion in the AWJSTOLR aircraft' 
‘simulator. 

6.2 CONCLUSIONS 

The quality of the flight data- provided was excellent. 
Inputs controlled by the avionics computer were exceptionally 
repeatable and accuracies associated with parameters were quite 
good for this reason. The data processing was nearly routine 
with the identification and run-to-run repeatability of the 
parameter estimates was good. 

The results indicate that the engine is represented, for 
outer loop simulation and control evaluation, as a second order 
system. Thrust outputs must be scheduled on fan and compressor 
speeds separated to model the forces generated by the engines. 
The model developed is accurate above about 18° throttle angle. 
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The behavior o£ the engine and control, aside from the hyster- 
esis, is quite linear. An insignificant asymmetry in the 
acceleration and rates was found for these powe.r levels. 

The engine-to-engine repeatability of the parameter 
estimates was within 10% except for the compressor time con- 
stant. The right engine time constant is 35% longer than the 
left. Fan speed responses are the same. The fuel flow to 
rotor speed gains are nearly the same, hence, the engine fuel 
consumptions are nearly equal. Fan speed responses are 50% 
slower than compressor time constants . 

The utilization of automated parameter identification pro- 
cedures has provided a significant improvement over classical 
methods. The number of parameters estimated is far greater 
than previously available from graphical techniques. Run- to- 
run repeatability has been demonstrated for many of the 
variables. 

The established system- identification routines can be used 
in the" design and evaluation of future flight test applications. 
Several of these areas” are mentioned below. 

The tests described in this report primarily produced 
linear dynamic results. The scope of the effort was restricted 
to generally linear models relating the engine response near 
full power. A more complete representation would use response 
data from a larger power envelope to estimate both steady state 
and dynamic coefficients. 

A nonlinear model can be used to match the response using 
this data for the full envelope of operation including nozzle 
angle and choke effects. ’Pressure data in the duct and tail- 
pipe could be matched much more accurately using a quasi- linear 
model of the following form: 

p = -X(p-fCN)) 
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•where the function, fC*)» describes the nonlinear static 
operating line response in pressure- The results for duct 
pressure response indicate that the choice of X is not im- 
portant since X > 20 sec ^ for both the duct and tailpipe 
dynamics . 

The accuracy of the identification procedure can be improved 
using optimal input design. These techniques build upon the 
results of this study. Estimates for model parameters developed 
from sawtooth and doublet inputs form the input design model. 
Throttle input time series are generated from the design pro- 
cedure which could be programmed on the AWJSTOLR aircraft flight 
computer as throttle rate commands. Data from these flight tests 
would be used with the identification program to calculate a 
revised set of parameter estimates. This procedure would produce 
far more efficient utilization of the flight test time by en- 
hancing parameter identif iability . Effects of nozzle position 
and choke deflection on thrust dynamics could be efficiently 
evaluated. 
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APPENDIX A 

LISTING OF SPEY ENGINE SIMULATION 
AND IDENTIFICATION PROGRAM 


PAGE BLANK NOT >XLMED ' 



97 j>r0CS2DIN<» page NOT 



pKlifaRAfl ! 

iHkV 76/7b Ol>T = l FTN 5+010 

10/06/77 

11.15 

1 


PRDGHAH SPEY (iNPUFfllUT put ♦TAPE5=INPUT*TAPE6=OU1PUTiTAP£2=/12QOO) 

BUM 7 

1 


C 


SPEY 

3 


C 


SPEY 

0 


c 

SPEY engine SIHULAUnN AND IDhNllPIC ATION PHOGRAtl 

SPEY 

5 

■i 

c 

AriJSTOLHA fLlUMT TEST PROGRAM 

SPEY 

6 


c 


SPEY 

7 


c 

systems control, INC (VT> 

SPEY • 

e 


c 

laoi-PAGE MILL ROAD 

SPEY 

9 


c 

PALO alto, CALIFORNIA 

SPEY 

10 

to 

c 

ATINI R.L.OEHOKf 

SPEY 

11 


c 


SPEY 

12 


c 


SPEY 

13 



external funcfunct 

BUN27 

1 



Real kth,kaf»kai,ka2,kfsf,kfijias 

SPEY 

15 

15 


LEVEL 2*X0ATA 

JUM5 

1 



common /engine/ 

SPEY 

16 



A pUINJ • TTHE » OELTC , DELTHS , RpHSN , 

SPEY 

17 



B XNHHSN • XNLMSN , PT6MSN , PTDM3N , 

SPEY 

le 



C DEUTPo i.wFHEST , XNH , XNL , PT6 ,• . 

SPEY 

19 

20 


0 PTO 1 OELTT , DELITS , hFMB , WFFLWM 

SPEY 

20 



E wFEST , OELTlt * dELTTE , EPSRPH , 

SPEY , 

21 



F DTHROT , DAREA , FHSl , FMS2 , FMFS1-, 

SPEY 

22 



6 FHFS2 , THRERR , pULERR , XH5ERR , 

SPEY 

23 



H XLSERR f PT66RR , PTDERR 

SPEY 

20 

25 


COMMON /engine/ OUMKIOOJ 

■JUNIO 

2 


c 


SPEY 

25 



common /CONSTS/ STAftT,CH(0)» 

SPEY 

26 



A OLINKliALO , ALl , ALE , RPOTO • 

SPEY 

27 

' 


tt PPOTl , DLINK2 , KTH , ALPHA , KAF , 

SPEY 

2B 

30 


C CPCHIR , ZETA , OH , OMN , SIG , OT , 

SPEY 

29 



0 CEAG , CLAOO , CLaGG , CF , 

SPEY 

30 



E AfL , CHL f flag , FLAGP , FLAGG , ' 

SPEY 

31 



F tlag , tlago , tlagg » tgain , 

SPEY 

32 



G OLAG , DLAGO , DLaGG , DgAIN , 

SPEY 

33 

35 


H KAi • XA2 , EnGNO , EPSBIS , 

SPEY 

30 



T KFSF , KFBIAS , Bit 

SPEY 

35 



COMMON /CONSTS/ 

SPEY 

*36 



A pTOBls , PTbBiS , XNLDI3 , 

SPEY 

37 



0 NM »NP, NF»N0P,NFIRST,NINC, SPRINT fCpRNT, 

SPEY 

33 

«0 


C lPARM(20>,ERRpNTt6) ,R(6 ),mSIG, 

SPEY 

39 



0 EPS , delta , MAXfN • lOPT » PARHMQ(O) * 

SPEY 

' 00 



E FILENO , STINE , FLGTNO , HUNDEL , HONOPO , 

SPEY 

01 



F hUNMFF , HUNNH , HUNNL r hUNPTO * HUNPTD 

SPEY 

02 



COMMON /CONSTS/ DUM2(100) 

JUNIO 

•3 

il5 

c 


SPEY 

03 



common /DATA/XOAIA(10,6000) • - 

JUL20 

1 


c 


SPEY 

05 



COMMON /SYM8OL/ND5yMB,NCSYHB,DSYMB(50) ,C3YM8(150) 

SPEY 

06 


c 


SPEY 

0 7 

50 


COMMON /PloTS/CPLOT*PLOTF, INOI, A (200,7) , imago ( 5 IBUiPDAT A (200,20) 

JUMS 

2 


c 


SPEY 

09 



DIMENSION P(20j , P ARM ( 1 ) , F (1 0 ) , X J AC ( 200 ) , X J T J ( 20 0 ) , WORK ( 250 ) 

SPEY 

50 



dimension HLAB(223 ,nElG(20) ,ZVECTt000),NK(250) 

JUL26B 

1 



dimension HESS(20) 

JUN26 

1 

55 


dimension 0ISP(200) 

JUN30 

1 


c 


SPEY 

52 



equivalence (PARH( I ) ,START) 

SPEY 

53 


PAlib I 



»'Kui;pam aptr 


76/76 ClHI = l 


{■TN a.s+ai'j 


10/06/77 11.15.06S 


I’AGt 


60 


65 


70 


75 


eo 


85 


<£> 

00 

90 


95 


loo 


Ids 


iTo 


lUHEUsT /INPUT/ StAIITiCH* 

A DLINKliA(,0»ALl»AL2tUPOTO»DPOTl«DUNK2tKTHiAl.PHAtKAFfCfCNTKi 

0 ZtTAiOHfOHN*SiGil)T»CLAGfCF> AFLtAHLiFLAGtTUAUtTGAlNiDLAG, 

C 0GA1N,KA| *KA2tENGMO*Ep38lSfPTpai3»PT6Bl3fXNL018f 

P NMfNP,NF»NDP*NFlR5l »Nl. AST. NINC iSPRINTf CPRNT t IPARHf ERRPNT ♦ 

E RiNSlr,.EP8iD£LTAiHAXFNiI0PT»PAKHHQiKFSF»KFBlA8fBIT» 

F FlLENO.STlHE.FLGlMO.HUNPfcL.HUNDPOftlUNMFFtHUNNHfHUNNLl 
6 HUNPT6»HUNPTP,PFAc.CPl.0T.10PnM,IEHESSt 
H HESS.cHt. 

C 

c 

data (RLADCI) rl=l»2)/'tHPARM»0HElER/ 

C 

C ♦♦♦INITIALIZE VAftIABLts^^+ 

C 

call setup 

C 

OFACaO, I 
CPLOTaO.D 
PLOTPaO.O 
lOPTINsO 
IEHESSsO 
C 

C ♦♦♦head input DATA+A+ 

C 

REAP(5, INPUT) 

CALL CPRINT 
C 

DT=0,02TN1NC 
C ♦♦♦read OATATAt 
C 

IF(BIT,LT,0.*0) call rdata 
IF(BIT,GE,0,0) call toata 
c 

C ♦♦♦iNSfcRT INITIAL PARAMETER GUESSES*** 

C 

DO 220 IstiNP 
INOalPARHfl) 

PI n=PARH(IND) 

220 CONTINUE 
C 

c ♦♦♦pfrtubb initial guesses for bit case*** 

C 

IFTBIT.LT.O.O) go to 200 
PO 225 Isl.NP 
P(n=P(I)*(l .APFAC) 

225 continue 

c 

200 CONTINUE 

IFtlOPTIM.EQ.U GO TO 235 
C 

c ♦♦♦optimization routine*** 

C 

CALL ZXSSO I FUNC*NF,NP( NS1G.EPS .pelt A. 

1 HAXFN.IoPT .PAHUHQiP.SSOiFi 

2 XJAC.NF .XJTJ.HiikK. INFFHiIER) 


6PEY 

SO 

SPEY 

55 

SPEY 

56 

SPEY 

57 

SPEY 

58 

SPEV 

59 

SPEY 

60 

SPEY 

61 

JUN28 

2 

AUGl 

1 

SPEY 

63 

SPEY 

64 

SPEY 

65 

SPEY 

66 

SPEY 

67 

SPEY 

60 

SPEY 

69 

SPEY 

70 

SPEY 

71 

SPEY 

72 

SPEY 

73 

JUN27 

3 

JUN2S 

4 

SPEY 

74 

SPEY 

75 

SPEY 

76 

SPEY 

77 

SPEY 

76 

SPEY 

79 

JUL28 

2 

SPEY 

60 

SPEY 

61 

SPEV 

62 

SPEY 

63 

SPEY 

64 

SPEY 

65 

SPEY 

66 

SPEY 

87 

SPEY 

86 

SPEY 

69 

SPEY 

90 

SPEY 

91 

SPEY 

92 

SPEY 

93 

SPEY 

94 

SPEY 

95 

SPEY 

96 

SPEY 

97 

.SPEY 

96 

SPEY 

99 

0UN27 

4 

SPEY 

100 

SPEY 

101 

SPEY 

102 

SPEY 

103 

SPEY 

104 

SPEY 

105 


^8 

sk 

Es^ 



H II t’ I ^ 


I..,/ 1 


I ■ 1 


^PW-l 


IN 


Oi<l| 


115 

fiO TO 230 

JUN27 

5 


C 

JUN27 

b 


c ♦♦♦QuaSI-NeWTON tlHT IMIZATION*** 

JUN27 
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IFtIER.EQ.132) MRlTE(6fll3) 


SPEY 

• 169 



IFtIER.EQ.133) MRlrE(6ill4) 


SPEY 

170 
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E FILENO » STIMt t FLGTNO , HUNDEU i HUNDPO , 
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F rONMFF t HUNHH 1 HUNNL • HUNPTS » HUNPTO 

ENGINE 

39 



c 


JUN29 

1 


i)0 


SOFT IXI»XMfXBJ=AMAXl C0,,XIf ((XI-XB)*X6)/(XB-XM)+XB) 

J0N29 

2 



c 


ENGINE 

40 




IF(8TART,EQ.0,0) Go TO 10 

ENGINE 

41 



c 


ENGINE 

42 



c 

♦♦♦INI IIALIZATION PASSt’M- 

ENGINE 

43 


as 

c 


ENGINE 

44 




THRERR = 0,0 

JUMOA 

1 




FULERR s 0, 

JUNHA 

2 




XKSERR 3 OiO 

JUNiaA 

3 




XLSERR s 0. 

JUNHA 

4 


so 


PT3ERH = 0,0 

JUNHA 

5 




PTOERR = 0,0 

JUNHA 

6 




1)LINKI = AL2*0ELTC**2FALI*DELR*AL0 

JOM7B 

1 




OL INK Is soft (01. INK 1,-1 00. »0, OS) 

JUN29B 

1 




oeltt=deltc-olinm 

JUM7B 

2 


55 


DELTISsOELTC 

ENGINE 

47 




DlHROTsOELll-OLlNKa 

JUM7B 

3 



c 


JUN29B 

2 



si3TTRWJ I !k^ 


“liiMir 




li 


Ul 


feO 


fcS 


70 


75 


BO 


O 

Ln 


85 


90 


95 


too 


105 


llO 


C ♦^stability limits*** 

.1UN29B 

3 

c 

JUN296 

U 

CLA(^=SoFT{CU^<J*-10o,,0.l ) 

JUN29U 

s 

FLAG=SOFT(FLA(5*-I0o.f6,l ) 
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50 
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CPRINT 

2A 



EQUIVALENCE (PARri(l)iSTART) 

CPRINt 

25 

25 


DATA (iNDKIJiialrlO)/ 

CPRINT 

26 



* lOf I6,ll,l2,13f lO,lS»21,22iA3/ 

CPRINT 

27 



data (iNDPUlf I=lf 10)/ 

CPRINT 

■ 28 



A l7»10,ie«20t06t07,08f22t2Af4a/ 

CPRINT 

29 



data (lND3Cn»l = l HO)/ 

CPRINT 

-30 

30 


♦ 26t32f 29,30f3l i52,27»28»33»30/ 

CPRINT 

31 



data (IND/|(I)|I = 1*10)/ 

CPRINT 

32 



♦ 55»30,51 ,39,il2tS0,36f 37f «0f 01/ 

CPRINT 

33 


c 


CPRINT 

30 



OAlA ENLA8/6HLEFT j6HRIGHT / 

CPRINT 

35 

35 

c 


CPRINT 

■36 



HR1TE(6i 100) 

CPRINT 

37 



100 FORHATdHli///) 

CPHINT 

36 


c 


CP.RINT 

39 



WRITE(6»101) 

CPHINT 

40 

/JO 


101 F0RMAT(50X»33HSPEY ENGINE SIMULATION PARAMETERS/ 

CPRINT 

41 



♦ 50X*25HAhJSTUL RESEARCH AIRCRAFT/ 

CPRINT 

42 



♦ 5AXi25HSYSTeHS CONTROL* INC (VT)/ 

CPRINT 

43 



♦ 57XT20H1801 PAGE MILL ROAp / 

CPRINT 

04 



♦ 56X*21HPALU ALTOt CALIFORNIA//! 

CPRINT 

45 

/<5 

c 

' ^ 

CPRINT 

46 



wr1TE(6« 102) 

CPRINT 

47 



102 F0HMAT(56Xf 21(1H*) ) 

CPRINT • 

46 


c 


CPRINl 

49 



IF(OIT.LT.O.O) GO TO 10 

CPRINT 

50 

50 

c 


CPRINT 

51 



HRITE(6*103) bir 

CPRINT 

52 



lOi FnRHAT(56Xf IH*. 19X, IH*/ 

CPRINl 

53 



♦ 56XtlH**3Xf UNBUILT IN TESTiSXilH*/ 

CPRINT 

.54 



»• 5f.X*lH*»3Xf5HRAMP=tF5.1f6H DEG / 

CPRINT 

55 

55 


56XflH*f 19X,1H+/ 

CPRINT 

56 



♦ 56Xi21(|HH) 

CPRINT 

57 



on TO 20 

CPHINT 

58 




bUHKUUIlWF. UFKlNI 76/7b UF1=J 
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t 

5 


Ifl 


15 


30 


35 



00 


05 


50 


subroutine OPHINI (M) • . CPH1N1 

CONNON /ENGINE/X(5«) DPRINT 

COHKON /SyU80i./N05YMB»NC3YHB»DsYMB(50) tCSYHBdSO) DPRJNT 

COMHON/CONST8/OOMtt50J DPHINT 

CflHNOU /PLOTS/cPLOT»PLOTF«lNDliA(300«7)»IHiGO(5lSU»PDATA(200*aO) JUM5 
C OPRINT 

0 DPRINT 

dimension PUOjSOJ CPRINT 

C DPRINI 

data InD/0/ OPRINT 

C DPRINT 

MDSYM0=35 DPRINT 

INDsIND+1 OPRINT 

DO 1 IstfNDSYNll OPRINT 

P(INDn)=xm DPRINT 

1 continue OPRINT 

C DPRINT 

!P(N,F.Q,0.ANI).IND«NE.10) OU to 200 CPRINT 

IFIPLOTF.RT.O.O) Go to 200 DPRINT 

C CPRINT 

WRITE(6»100) DPRINT 

too FORMATtlHl i09X*22HSPEY ENGINE SIMULATION*/ DPRINT 

♦ 5ox,3'jhengine Parameter identification proghahi/ cprint 

♦ 50X»20HSVSTEHS C0NTR0L»INC (VT)i/ DPHINT 

♦ 50X,20HpALO ALTO»CALIFORNIa*/ dprinT 

♦ sox, 1 INVERSION OOli/) OPRINT 

C DPRINT 

DO 2 IsIfNDSYHU DPRINT 

RRllEIbflOl) n5YHH(I)»(P(j,n»J®l»IND) ' OPRINT 

a continue DPRINT 

lOl ro«MATClXtA'6f 1X»10(6T2,6)J DPRlNt 

IND=0 OPRINT 

C ^ ' DPRINT 

c ♦♦♦Save plot data if plotting»*+ dprxnt 

c DPRINT 

200 continue DPHINT 

IFIPLOTF.EO.D.O) return DPRINT 

C DPRINT 

INOlBluDl+1 DPRINT 

IFdNOl.GT.aOO) RETURN DPRINT 

C , , OPRINT 

P0ATA(lNDdl)=X(2) OPRINT 

DO 210 I=h7 DPRINT 

PDATA(lNDli2^n=X(i+2) OPRINT 

2 to continue cprint 

C CPRINT 

cn 220 1=1*7 OPRINT 

PDATAdNDlf 2'l‘H-n=X(9+n cprint 

220 continue DPRINT 

C ' OPRINT 

return OPRINT 

END . CPRINT 
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C 


CPRINT 

59 



10 CONTIMIIE 

. CPRINT 

60 

<<0 

C 


CPHINT 

61 



IbENGMO 

CPRINT 

62 



JsFLGTNO 

CPRINT 

63 



KsFIUEnU 

CPHINT 

6A 



LbSTIHE 

CPRINT 

65 

65 


WRlTE(6tlO^> ENLAB(I)*JfK,L 

CPRINT 

66 



roa F0HKAT(56X* IH*» 19X, lit*/ 

CPRINT 

67 



♦ S6XilH*,2X»A5t ;H ENGlNF.»i|Xf2H */ 

JUN17A 

A 



♦ 56Xiih*iih i15h Flight number lUfiH*/ 

JUM7A 

5 



* 56X( 1H*(I>X« 1 IHI^UN number fI2t3Xf2H */ 

0UM7A 

6 

70 


♦ S6XflH*.2XijHSri*113f2H ♦/ 

JUM7A 

7 



♦ 56X»2ltlM*)) 

CPRINT . 

72 


C 

. 

CPRINT 

73 



20 CONTINUE 

CPRINT 

7A 


C 


CPHINT 

75 

75 


, W«lTE(6»t051 

CPHINT 

76 



105 FORMAT(26Xtai(lH*3/26Xt 1H*» 19X, lH*f 19X»IH*»19X» lH*il9Xf IH*/ 

CPHINT 

77 



* ?6Xf lH*f5X»aHTHR0TTLE|6Xf lH»»7Xf7HCONTROL. 

CPRINT 

78 



* 5Xf IH*»7 Xi6HENGINe» 6X, lH*,i|Xi 

CPRINt 

79 



* l2HDUcT/EXHAUsTt5XilH*»/ 

CPRINT 

60 

60 


* r 6X, lH*i l9X,lH*t l9X»lH**iqXtlH*fl9XilH*/26Xf 8K1H*)) 

CPRINT 

81 


C 


CPRINT 

62 



on 50 lalilO 

CPRINT 

63 



n = iNDi(n 

CPRINT 

6A 



r2=lND2U) ■ 

CPRINT 

65 

65 


I3»IN03( I) 

CPRINT 

66 



ia=iND6a) 

CPRINT 

87 


C 


CPRINt 

68 



URU£(6»1103 C3VMB(U) ,PARM(in* 

CPRINT 

89 



♦ CSYN«(I2) iPARM(I 2) f 

CPRINT 

90 

90 


♦ CSYMB(IJ!)iPARH<I5)* 

CPRINT 

91 



♦ CSYI1Q(I«) iPARHtm) 

CPRINT 

92 


c 


CPRINT 

93 



110 format (26Xf 6 ClH*fA6ilH=»G12,6) , 1H*J 

CPRINT 

9A 



50 COimNUF. 

CPRINt 

95 

95 

c 


CPRINT 

96 



WRnE(6»llU 

CPHINT 

97 



111 F 0 RMATr 2 t>XteitlH'O//J 1 

CPHINT 

98 



RRlTE(6il20) 

CPRINT 

99 



120 FnRMAT(26X*9l(im‘3 (bXiiadH*)) 

CPRINT 

100 

loo 


URlTE{hfl2t) 

CPRINT 

101 



121 F0RMAT(?6Xf iHd 12X, 16HFULL POWER POINT t 1 1 X • IR* > 6Xi IH* • AX • 

CPRINt 

102 



* SiHOPTIHiZATION PAHAHETEPSfSXf IH*) 

CPRINT 

103 


c 


CPRINT 

lOA 



WRITE(6* 1221 

CPRINT 

105 

105 


122 FORHAT(26X*lH*f39X,lR*,6Xf lHAii2X»lH*l 

CPHINT 

106 


c 


CPRINT 

107 



WRITEIei 131 1 HUNDE|,»NM 

CPRINT 

106 



UR1TE(6i 132) HllNOPOfNP 

CPRINT 

109 



wtJlT£(6f 1331 HUNMFFiNF 

CPRINt 

110 

110 


WR1TE(6i IX fll hunnh.ndp 

CPRINT 

ill 



HRn£t6*l35) HUNNL,NF1RS1 

CPRINT 

112 



WRlTE(b*l36) HUNPT6 »NInC 

CPRINT 

113 



IPRNT=SPRINT 

CPHINT 

llA 



wRlthlftf 1 573 HUNPTOiIPRNT 

CPHINT 

115 



116 


SUB‘<UUfINE 1 

LPiUi'JI 76/7b OPT^l FTN 

4.5+414 
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11.15 

tl5 

C 

• 


CPRINT 

1 16 



IH FORMAT! 


CPfilNT 

117 



♦ 26XtlH’F»lXM9HTHR0TTLE POSITION s,g10,0*9H 

DEGREES »IH*» 

CPRINT 

lie 



♦ hX»lH*t2?H' NO,OF HtASUPEMENTS » tIOfSXilH*) 


CPRINT 

119 



112 FORMAT! 


CPRINT 

120 

120 


♦ 26X» IH*»1X* lOHTHROTTl-E PICKOFF S|G10,0|OH 

COUNTS tlHt, 

CPRINT 

121 



♦ bX.lH+taiH NO, OF FAHAMtTEBS c tIOiSXilH*) 


CPRINT 

122 



113 FORMAT! 


CPRINT • 

123 



♦ abX^lHtilXMRHHAIN FUEL FLOW =»G10,0#<>H 

LBM/HOUR»lH*i 

CPRINT 

124 



♦ ftX,lHt.23H MO, OF error TERMS = iIOtSXtlH*) 


CPRINT 

125 

125 


llA FORMAT! 


CPRINT 

126 



* 26X* IH+i IX* I’JHCOMPRESSOB SPEED a*Gl0,0|9H 

REV/MIN ilH'+i 

CPRINT 

127 



♦ >X*IH*.23K NO, OF data POINTS = iHltSXflH^l 


CPRINT 

126 



115 FORMAT! 


CPRINT 

129 



♦ 26X*1H», IX.19HFAN SPEED SfGlO.OiRH 

REV/HIN tlH*» 

CPRINT 

13Q 

130 


♦ AX,1H»*23H NO, OF FIRST POINT c ,I0,SX|1H*) 


CPRINT 

131 



Ufa FORMAT! ■ 


CPRINT 

132 



♦ 26X*lH»t!Xfi9HAUGri£NT0R PRESSURES. G10,«*9H 

LBF/SQINf IH** 

CPRINT 

133 



♦ 6XflH*i23H NO, OF INCREMENT * tlA.SXtlH*) 


CPRINT 

. 134 



U7 FORMAT! 


CPRINT 

135 

135 


♦ • 26X.1H*, 1X*19HDUCT PRESSURE s,G10,i|*9H 

LBF/SQIN*IH*| 

CPRINT 

136 



♦ 6X,lH*f25H DATA PRINT INTERVAL’ tIi||5X*lH+J 


CPRINT 

137 



WRlTEISiHO) 


CPRINT 

136 



lOa F0RMAT!26X.1H+,39X,1H*»6X, jH*i32Xf IH*) 


CPRINT 

139 


C 



CPRINT 

140 

UO 


ICPRNTsCPRNT 


CPRINT 

141 



WRlTE!6*t«lHCPRNT 


CPRINT 

142 



lai FORMAT(2faXf01!lHt),6XflH*»2lK CONSTANT PRINT FLAG: 

>.2XiI4*5XilH») 

CPRINT 

U3 


C 



CPRINT 

144 



WRITE(6.1021 F.PS 


CPRINT 

145 

105 


WRITE(6*103) delta 


CPRINT 

146 



w«lTt!6*lA9) Maxfn ; 


CPRINT 

147 



WRITE(S«105] NSIG 


CPRINT 

146 



U2 FOHMATI73X* IH*, 


CPRINT 

149 



♦ 2lH epsilon s*GlI,5ilH>t') 


CPRINT 

150 

150 


1113 FDRHATlTSXtlMtf 


CPRINT 

151 



♦ 2lH delta sjGlI.S.lH’F) 


CPRINt 

152 



F0RHAT!73X*1M** 


CPRINT 

153 



* 21H MAX FUNCTION EVALS a t 2X » 14 i 5X * IHT) 


CPRINT 

154 



US FORHAT(73X* lH+» 


CPRINT 

155 

155 


♦ 2JH NO, OF SIGnIF DIGITSs * 2X f 14 *SX * IH*J 


CPRINT 

156 


C 



CPRINT 

1S7 



WRlTE!faiI50) 


CPRINt 

ise 



ISO FORMAT! 73X.3R! INI')) 


CPRINT 

159 


C 



CPRINT 

160 

l60 


RETURN 


CPRINT 

161 



end 


CPRINT 

162 


HAGt S 
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!)A^ 




7 ^> 


^ ^ 


10 


11 


5,0l 


1 


BLOCK data SYHHBLS - : 

symbols 

2 



CoUHON /8YMUOL/NO8YHBiNCSVMB»OsyH8(50) tCSYHBClBO) 

SYMBOLS 

3 


C 


SYMBOLS 

H 



DATA NpSYMD»NCSYHtl/35i 150/ 

SYMBOLS 

5 

5 


data OSYHBtCSYlta/ScOl'OH / 

SYMBOLS 

6 


C 


SYMBOLS 

7 



data (DSYHB(n»I = l ilO)/ 

SYMBOLS 

6 



A 6HP0INT , 

SYMBOLS 

5 



B 6 H 1 H 1 C 1 

SYMBOLS 

10 

10 


C 6 HDEUTC t 

' SYMBOLS 

11 



0 *hdelths» 

SYMBOLS 

12 



E AHWFHSN , 

SYMBOLS 

13 



F 6 HXNHMSN 1 

SYMBOLS 

H 



G 6 MXNLMSN 1 

SYMBOLS 

15 

15 


H 6HPT6mSN» 

SYMBOLS 

16 



I hHPTDMSNt 

SYMBOLS 

17 



J 6HDELTP0/ 

SYMBOLS 

18 



data (DSYMB(I)iI=il»20)/ 

SYMBOLS 

15 



A 6HWF.MEST. 

SYMBOLS 

20 

20 


B 6i<XMH ' f 

SYMBOLS 

21 



C 611XNL t 

SYMBOLS 

22 



0 6HPT6 , 

SYMBOLS 

23 



E 6HPTO ♦ 

SYMBOLS 

21) 



F 6H0ELTT 1 

SYMBOLS 

2S 

25 


G 6HDELTTS. 

SYMBOLS 

26 



H 6HWFHR , 

SYMBOLS 

27. 



.1 6KHFFLWH» 

SYMBOLS 

28 



J 6HWFEST / 

symbols 

29 



DATA (DSYMBd) »I“2i r30)/ 

SYMBOLS 

30 

30 


♦ 6BDEUTT1* 

' SYMBOLS 

31 



A 6H0tLTT2, 

SYMBOLS 

32 



B 6HEP3rPHi 

SYMBOLS 

■33 



C 6 HDTHROT 1 

SYMBOLS 

34 



0 6HDAREA * 

SYMBOLS 

35 

35 


E 6HFHS1 , 

SYMBOLS 

36 



F 6HFMS2 ♦ 

SYMBOLS 

37 



G 6HFMF31 • ■ ■ 

SYMBOLS 

38 



H 6HFBFS2 1 

SYMBOLS 

39 



I 6HTHREBR/ 

SYMBOLS 

40 

<10 


data (DSYHafl)*I = 3l»<lO)/ 

SYMBOLS 

91 



. ♦ 6HFULeBR» 

SYMBOLS 

42 



A 6HXH8ERRt . 

SYMBOLS 

43 



a ohxlserri 

SYMBOLS 

44 



C 6HPT6ERR* 

SYMBOLS 

45 

05 


0 f.hptderr» 

SYMBOLS 

46 



E 6H , 

SYMBOLS 

47 



F 6H f 

symbols’ 

48 



G fcH , 

SYMBOLS 

49 



H bH , 

SYMBOLS 

50 

50 


t 6H / 

SYMBOLS 

51 



DATA (DSYMB(IJfI=4|»50)/ 

SYMBOLS 

52 



A 6H , 

SYMBOLS 

S3 



B ftH , 

SYMBOLS 

54 



C 6H • 

SYMBOLS 

55 

55 


0 6H , 

SYMBOLS 

56 



E , 

SYMBOLS 

57 



F 6H , 

SYMBOLS 

58 
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G 6H , 

SYMBOLS 

59 


H 6H • , 

SYMBOLS 

60 

bO 

I 6H , 

SYMBOLS 

61 


J 6H / 

symbols 

62 


data (c8YMB(1)i1=1 *10J/ 

SYMBOLS 

63 


A 6HSTART • 

SYMBOLS 

64 


B bHCHll) , 

SYMBOLS 

65 

65 

C 6HCH(2) I 

SYMBOLS 

66 


0 6HCH(3) « 

SYMBOLS 

67 


E 6HCH(4J • 

SYMBOLS 

6B 


F 6HCH(5) f 

SYMBOLS 

69 


G 6HCHt6) 1 

SYMBOLS 

70 

70 

H 6HCHC7) t 

SYMBOLS 

71 


I 6HCH(8) . 

SYMBOLS 

72 


J 6H0LINKI/ 

SYMBOLS 

73 


DATA (cSyMB(I)tUlltaO)/ 

SYMBOLS 

74 


A 6HAL0 « 

SYMBOLS 

75 

75 

a 6HAL1 , 

SYMBOLS 

76 


C 6HALH , 

SYMBOLS 

77 


D 6HBPOTO » 

SYMBOLS 

76 


E 6HDPOTI . 

SYMBOLS 

79 


F 6HPLINK2i 

SYMBOLS 

60 

no 

G 6HKTK • 

SYMBOLS 

61 


H 6HALPHA , 

SYMBOLS 

. 82 


I 6HKAF , 

SYMBOLS 

63 


J bHCFCNTR/ 

SYMBOLS 

64 


DATA <cSYHBmils2j*30)/ 

SYMBOLS 

65 

S5 

A 6H2ETA t 

SYMBOLS 

86 


B 6H0H , 

SYMBOLS 

67 


C bHOtlN , 

SYMBOLS 

68 


0 6HS1G , 

SYMBOLS 

69 


E 6HDT , 

SYMBOLS 

90 

90 

F 6HCLAG ♦ 

SYMBOLS 

91 


G 6HCLA6D t 

SYMBOLS 

92 


H 6HCLARC , 

SYMBOLS 

93 


1 6HCF , 

SYMBOLS 

94 


J 6HAFL / 

SYMBOLS 

95 

95 

DATA (CSYH0n)*I = 3l»«O)/ 

SYMBOLS 

96 


A 6HCRU , 

SYMBOLS 

97 


B 6HFUAG , 

SYMBOLS 

98 


C 6HFLAGD * 

SYMBOLS 

9.9 


0 6HFLACG ♦ 

SYMBOLS 

100 

loo 

E 6HTLAG , 

symbols 

101 


f 6HTLAgD t 

SYMBOLS 

102 


G 6HTLAGG , 

SYMBOLS 

103 


H 6HTGAIN ♦ 

SYMBOLS 

104 


I 6HDLAG , 

SYMBOLS 

105 

105 

J 6H0LAGD / 

symbols 

106 


OA7A tcSYHQm«l = 'll*5l)/ 

SYMBOLS 

107 


♦ 6H0LAGG , 

SYMBOLS 

108 


A 6H0GAXN , 

SYMBOLS 

109 


B 6HKAI- , 

symbols 

110 

no 

C 6HKA2 , 

SYMBOLS 

111 


0 6HENGNO t 

symbols 

112 


E 6HEF$nl5t 

SYMBOLS 

113 


* bUKFSj- , 

SYMBOLS 

114 


♦ hHKFHjAS, 

symbols 

115 



119 


iSTTiT orrsT 


“Till 4) 4 


t<nrnr daia uynuu 


t (JV li'fc / 7 / "TTTlb.flT^ "TSTSfc 


US 


120 


125 


130 


135 


lilO 


M5 


150 


155 


IhO 


165 


170 


♦ 6HBIT , 

SYHB0L3 

116 

f 6HPTO0IS, 

SYMBOUa 

m 

G,fcHPT6»IS/' 

SYHBOLS 

lie 

riATA (csynuaif l = S2f61)/ 

SYP80L8 

119 

H 6HX6iLflISt 

SYKliOLS 

120 

I , 

5YHBOL3 

121 

J 6HNP t 

SYFPOLS 

122 

A 6HNF . 

3YF601.S 

123 

0 6HHUH , 

SYPCOLS 

12A 

C 6MNFlRSTt 

SYFBOUS 

125 

D 6MNIHC 1 

SYMBOLS 

126 

E iHSPfilNTf 

SYMBOLS 

127 

F'6HCPRfJT f 

SYPQOLS 

126 

G 6HPARI11 f 

symbols 

129 

OAfA tCSYMBdJ f I = 62r71)/ 

SYMBOLS 

130 

H 6HPAHH2 t 

SYMBOLS 

131 

I 6HPAKH3 • 

SYMBOLS 

132 

J 6HPARm« ♦ 

symbols 

133 

A 6HPARH5 . 

Symbols 

134 

B <,tiPARH6 . 

symbols 

135 

C 6HPARH7 » 

symbols 

136 

0 6HPARM6 « 

symbols 

137 

E 6HPARH9 f 

symbols 

138 

F 6HPARH10, 

SYMBOLS 

139 

G tHPARHll/ 

symbols 

140 

data (cSYMBd) I l = 72«8t)/ 

SYMBOLS 

141 

H 6HPARH12. 

SYMBOLS 

142 

I 6HPARH13* 

SYMBOLS 

143 

J frHPARMl4i 

SYMBOLS 

144 

A 6HPARh15i 

SYMBOLS 

145 

8 <,HPARM16. 

symbols 

H6 

C 6HPARM17* 

symbols 

147 

0 68PAR»tl8i 

symbols 

148 

E 6HPAR»19» 

symbols 

149 

F 6HPARn20f 

SYMBOLS 

150 

G 6HERR1 / 

SYMBOLS 

1st 

DATA (cSYMBdJ |I = U2»9l)/ 

SYMBOLS 

152 

H 6HERH2 , 

SYMBOLS 

153 

I 6HERR3 f 

SYMBOLS 

154 

J 6HEHRq , 

symbols 

155 

A 6HERR5 1 

SYMBOLS 

156 

B 6HEHR6 • 

SYMBOLS 

157 

C 6HFACT1 , 

SYMBOLS 

158 

1) 68FACT2 t 

■ SYMBOLS 

159 

E. 6HFACT3 t 

SYMBOLS 

160 

F 6HFACT0 , 

SYMBOLS, 

161 

G 6HFACT5 / 

SYMBOLS 

162 

DATA (CSYMfld)|I=92f 101)/ 

SYMBOLS 

163 

M 6HFACT6 t 

SYMBOLS 

164 

1 6HNS7G f 

SYMBOLS 

165 

J 6HEPS , 

SYMBOLS 

166 

A 6HDEI.TA , 

SYMBOLS 

167 

b 6 H»^AXfM t 

SYMBOLS 

16B 

C 6HIOPT , 

symbols 

169 

D 6HPARMQ1, 

SYMBOLS 

170 

E 6HPARH02. 

SYMBOLS 

171 

F hHPARM03» 

SYMBOLS 

172 
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' HAtiE 

1 

subroutine ROATA 

■ 

JUM5 

6 



C 


JUMS 

7 



C ♦♦♦subroutine to RF.AO flight DATA**t . 

JUN15 

B 



c 


JURIS 

0 


5 

LEVEL ZfXDATA 


JURIS 

10 



COMMON /OATA/XoaTA(IO» 1) 

\ 

JUMs 

11 



C 


JUM5 

12 



common /CONSTS/ STARTiCH(6)i 


JUMS 

13 



A oLINKIjALO t AH t AL2 f BPOIO f 

JUMS 

H 


10 
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APPENDIX B . 
FLIGHT TEST PLAN 


JL, FI 7 towards Crows Landing and climb to 7500 ft. Estab- 
lish communication and telemetry with NASA2 as soon as pos- 
sible. If unable to establish communication, start the 
throttle sweeps when 15 DME west of Crows Landing. 

Z. Establish the airplane in the STOL mode on a -1-1/2° 
fit path angle. Once stabilized at a speed of LAS + 5 hts 
proceed with the following matrix of throttle rate square 
waves . 


-TT-T 

SQUARE 

WAVE 

DURATION. 

AMPLITUDE . 

RECORD LENGTH . 


— 

EH 

5 sec 

1.25° 

1/2°7sec 

• 10 

30 sec- 

II 

2.5^ 

l°/sec 

20 

$1 

11 

5.0° 

2°/sec 

40 

n 

II 

7.5° 

3°/sec 

60 

ti 


4 to 5 min + climb 

3, Climb back to 7500 ft. and proceed through the following 
matrix. 


DURATION 

GNJ 

AMPLITUDE 

GNI 

RECORD LENGTH 

10 sec 

100 

1/2°/sec 

10 

40 sec 

w 

11 

1.5°sec 

30 

11 


3 min + climb 

4. Climb back to 7500 ft and do the following case: 


DURATION 

GNJ 

AMPLITUDE 

GNI 

RECORD LENGTH 

ll^liyillllllll 

500 

.25°/sec 

5 

2 cycles 


2 min + climb 
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5. 


Climb back to 7500 ft. and do the following doublets. 


DURATION 

AMPLITUDE 

GNJ 

GNI 

RECORD LENGTH 

10 sec 

o 

o 

o 

20 

30 sec 

11 

" 1-1/2° 

30 

II 

II 

" 2° 

40 

II 

II 

" ■ 4° 

80 

II 

II 

" 7-1/2° 

150 

II 
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TIME 

Figure C 
Figure C 
Figure C 
Figure C 
Figure C 
Figure C 
Figure C 

Figure C 

Figure C 

Figure C 

Figure C 


APPENDIX C 

HISTORY MATCHES FOR PARAMETER IDENTIFICATION RUNS 

. 1 Throttle Linkage Match - Record 1, Left Engine 

.2 Throttle Linkage Match - Record 5, Left Engine 

, t' 

.3 Throttle Linkage Match - Record 1, Right. Engine 

.4 Throttle Linkage Match - Record 5, Right Engine. 

.5 Tailpipe/Duct Pressure Match - Record 1, Left Engine 

.6 Tailpipe/Duct Pressure Match - Record 1, Right Engine 

.7 Engine Parameter Estimation Match - Record 3, Right 
Engine 

.8 Engine Parameter Estimation Match - Record 7, Right 
Engine 

.9, Engine Parameter Estimation Match - Record 7, Left 
Engine 

.10 Fuel Control Estimation Match - Record 3, Right 
Engine 

.11 Fuel Control’ Estimation Match - Record 7, Right 
Engine 
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Figure C.7 (Continued) 
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Figure G.9 (Continued) 
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Figure C.IO Fuel Control Estimation Match - Record 3, Right 
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Figure C.ll Fuel Control Estimation Match - Record 7, Right 
Engine 
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